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ABSTRACT 

Experimental  results  for  the  complex  dielectric  constant  of 
the  air-cellulose-water  mixture  as  found  in  Douglas  Fir  and  Western 
Hemlock  are  given  as  a  function  of  moisture  content  and  temperature 
at  2450  MHz  -with  a  maximum  error  of  +  6$.  Theoretically  the  diel¬ 
ectric  constant  of  wood  can  be  predicted  to  +  10$  by  assuming  a 
Gaussian  distribution  of  relaxation  times.  Experimental  methods 
are  developed  and  explained  in  detail.  Present  dielectric  mixture 
relations  are  briefly  reviewed  and  a  more  general  case  of  multi¬ 
phase  mixtures  with  confocal  ellipsoidal  shell  inclusions  is  ana¬ 
lyzed.  Values  for  the  complex  dielectric  constant  of  bound  water 
as  found  in  cellulose  mixtures  are  estimated  on  a  semi-empirical 
basis  and  are  found  to  vary  with  moisture  content. 
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CHAPTER  1 


1 


INTRODUCTION 

Multiphase  dielectric  theory  is  important  in  many  aspects  of 
material  research.  One  area  where  material  research  is  extremely 
important  is  in  the  application  of  microwaves  as  a  source  of  highly 
refined  heat  energy  and  in  measurement  techniques  used  in  microwave 
processing  of  foods ,  chemicals,  plastics,  etc,,  and  microwave  mois¬ 
ture,  thickness  and  fault-detection  measurements.  To  do  these  jobs 
effectively,  efficiently  and  economically,  the  microwave  properties 
of  the  materials  used  must  be  known.  This  thesis  extends  the  dielec¬ 
tric  theory  of  multiphase  mixtures  with  application  to  the  wide  class 
of  cellulose-air-water  mixtures  as  found  for  example  in  wood. 

Given  a  dielectric  mixture  made  up  of  different  components  each 
with  its  own  dielectric  characteristics,  can  the  dielectric  behavior 
of  the  mixture  as  a  whole  be  predicted.  To  answrer  this  problem,  die¬ 
lectric  mixture  theory  was  developed  as  far  back  as  the  late  19th  cen- 

(29) 

tury  when  Maxwell  (1891)  tackled  the  problem  of  obtaining  the 
electromagnetic  characteristics  of  simple  dielectric  mixtures.  Since 
then  many  scientists  have  modified,  added  to,  and  improved  on  the 
early  theories. 

An  exact  solution  of  the  problem  is  possible  for  only  a  few 
geometrically  Y;ell-defined  systems.  For  many  naturally  occuring 
dielectrics  the  geometries  of  the  mixtures  and  the  interaction 
between  the  inclusions  can  only  be  described  statistically  at  best. 


' 


\ 


2 

In  the  case  of  cellulose-sir-water  systems,  as  studied  in  this  thesis, 
even  a  statistical  description  is  extremely  difficult  to  obtain  and 
only  a  relatively  crude  mathematical  model  of  the  macroscopic  mixture 
can  be  set  dovm  which  nevertheless  is  shown  to  be  a  useful  analytical 

tool* 

( 13 ) 

Development  of  dielectric  mixture  theory  has  shown v  '  that  par¬ 
ticle  interactions  and  statistical  variations  in  the  mixture  have  only 

a  second  -  or  higher  -  order  effect  and  that  the  particle  geometry  in 

(14) 

the  mixture  is  a  far  more  important  parameter  , 


as  a  result, , 

existing  averaging  laws  are  used  in  this  work  but  a  more  general 
inclusion  geometry  is  assumed  namely  that  of  ellipsoidal  shells  which 
leads  to  a  new  mixture  relation* 


The  problem  of  finding  the  internal  to  external  electric  field 
ratios  for  ellipsoidal  shell  inclusions  is  solved  here  for  the  first 
time.  The  final  solution  simplifies  in  special  cases  to  that  of 
ellipsoidal  and  spherical  inclusions.  Chapter  2  is  devoted  to  the 
development  of  the  new  mixture  relation,  A  brief  history  of  mixture 
theories  is  presented  arid  summarized  in  a  tabular  form  which  can  be 
used  as  an  easy  reference  chart. 


Many  mixed  dielectrics  exhibit  a  complex  functional  dependence 
on  temperature  and  moisture  content.  This  nonlinear  functional  be¬ 
havior  is  evident  when  trying  to  measure  the  moisture  content  in 
materials  using  microwave,  capacitance  or-  resistance  type  moisture 


-X 


. 

* 


meters  .*,f  The  interaction  of  temperature  and  frequency  dependence 
of  the  material  properties  makes  it  very  difficult  to  calibrate  these 
instruments.  To  analyze  these  functional  variations s  chapter  3 
develops  the  theory  of  the  temperature  and  frequency  behavior  of  those 
dielectrics  which  exhibit  a  distribution  of  relaxation  times.  The 
theory  gives  an  analytical  description  of  the  variation  of  dielectric 
constant  with  temperature  frequency  and  moisture  content  and  there¬ 
fore  is  a  useful  tool  in  the  analysis  of  those  systems  in  which  the 
dielectric  properties  are  crucial  in  determining  systems  efficiency ^ 
stability  and  economics. 

To  test  the  theories  presented  it  ?/as  necessary  to  make  a  large 
number  of  measurements  on  an  actual  complex  dielectric  mixture.  Two 
wood  species  Douglas  Fir  and  Western  Hemlock  were  chosen  because 
of  their  economic  importance  on  this  continent  and  because  wood  is  a 
good  example  of  a  much  v/ider  class  of  cellulose-air-water  mixtures. 
Chapter  4  describes  in  detail  the  microwave  measurement  method  chosen 
and  an  error  analysis  needed  to  obtain  percentage  errors  is  developed 

V 

in  Appendix  A.  Overall  measurement  errors  are  shown  to  lie  between 
+  5  and  +  10#  eventhough  the  basic  technique  used  is  accurate  to 
better  than  +  2#. 

Results  of  the  measurements  on  Douglas  Fir  and  Western  Hemlock 


-x-  These  meters  are  commonly  used  in  the  forest  product  industry  and 
their  reliability  and  accuracy  has  been  much  debated  in  the  technical 
world . 


1%.  * 
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are  given  in  chapter  5  and  Appendix  B.  It  is  shewn  that  the  temperature 
and  frequency  behavior  of  the  dielectric  properties  of  wood  at  varying 
moisture  contents  can  adequately  be  described  by  the  theory  given  in 
chapter  5.  Agreement  between  the  experimental  and  semi-theoretic al 
behavior  of  wood  is  better  than  +  10$, 


Chapter  6  is  devoted  to  the  application  of  the  new  mixture  re¬ 
lation,  developed  in  chapter  2,  to  the  wood  mixture.  Since  the  ad¬ 
sorbed  water  in  wood  is  molecularlv  bound  in  varying  degrees  to  the 
basic  wood  structure  the  dielectric  properties  of  this  adsorbed  water 
differ  considerably  from  the  pure  liquid  water  properties.  Hence, 
the  new  mixture  relation  was  used  to  s emi -empiric ally  obtain  estimates 
of  the  dielectric  properties  of  the  bound  water  as  a  function  of 
moisture  content,  data  previously  not  available.  The  model  assumed 
for  the  wood  mixture  below  fiber  saturation  is  one  containing 
ellipsoidal,  needle-like  air  cavities  and  thin  shell-type  water 
inclusions  dispersed  in  a  homogeneous  cell-wall  material.  Calcu¬ 
lations  based  on  this  model  agree  favorably  with  experimental  data 
on  Douglas  Fir  and  Vie  stern  Hemlock. 


w"  X 
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CHAPTER  2 
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DIELECTRIC  MIXTURE  THEORY 

I 

Introduction 

Relating  the  dielectric  constant  of  a  mixture  satisfactorily  to  the 
dielectric  constant  of  its  component  has  been  a  major  problem  in  dielectric 
theory.  The  complete  solution  awaits  mainly  the  development  of  an  accurate 
theory  of  elementary  particle  interaction  based  on  statistical  mechanics 
which  is  by  no  means  a  trivial  task.  A  brief  review  will  be  given  here 
indicating  the  basic  mixture  problem  and  some  of  the  contributions  made  in 
the  literature  towards  its  solution.  A  generalized  approach  will  then  be 
given  to  the  problem  of  mixtures  following  the  work  of  Taylor One  of 
the  more  important  parts  of  any  mixture  theory  is  to  find  the  relationship 
between  a  particle's  internal  electric  field  and  the  average  field  present 
in  the  remainder  of  the  mixed  dielectric.  A  general  solution  to  the  elec¬ 
trostatic  problem  of  an  ellipsoidal  shell  will  be  given  with  a  view  to  the 
solution  of  the  macroscopic  cellulose-air-water  system  to  which  the  diel¬ 
ectric  mixture  theory  will  be  applied. 

A  chart  of  the  important  mixture  relations  is  presented  in  Table  2.1. 

During  the  course  of  this  research  an  excellent  review  of  mixture  relations 

(64) 

was  published  by  Van  Beek  ,  summarized  in  tabular  form.  However,  Table 
2.1  given  here  presents  the  major  mixture  relations  from  which  all  the 
others  can  be  derived.  It  further  covers  the  work  of  Taylor^  and 
Franck ^ 61  ^  both  of  whom  have  made  significant  contributions  by  stating 
the  mixture  problem  and  solution  in  more  general  terms. 


I 


fc:  1 


Table  2.1  -  Summary  of  most  important  macroscopic  mixture 
relations 
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■Inclusion 

shape 

Mixture  Formula 

Conditions  and 
References 

Ellipsoids 

2-phase 

£*=  eh 

All  interaction  effects 
assumed  negligible; 
vlaxvvell  ( 29  j  Fricke. 

Arbitrary 

2-phase 

(£m-£\n)/(£n 1+W)=  Vi  (El  + 

- ....  - 

Semi-empiric  al; 

Wiener*-^) . 

Arbitrary 

2-phase 

(£„-£h)ferQ  =  £'=L(i-'')''il/Hv‘) 

k  is  an  empirical  factor 

- - , 

Empirical;  Pearce  (^3)^ 

Spheres 

2-phase 

£m  -  £h  =  5V;  £h /f[(£P*2£h)/(£r£j>V;}  . 

Interactions  allowed 
for  by  ;  Leivin(^)  . 

Needles 

Discs 

Spheres 

2-phase 

'-Vi  =  {(£i-£„)/£;-£lj}.(2£;  +£h)/(2£:  +£w)  ; 

i-Vjr  {(£;-£w)/ft-4 ))• 

Interactions  partially 
accounted  for  in  the 
ierivation;  . 

3ruggeman-Niesel'‘-L‘i'' . 

• 

Arbitrary 

multiphase 

M  1 

bi£M  =  V,  L«e,n  +5.  ViLn£i  Empirical; 

lat|  ;Lichtenecker(14,21) . 

i  _ 

Spheres 

2-phase 

{ (W  £h) /3£m)  =  V; (£;-4i/(£;  +2.£m) 

Short  range  statistical j 
variations  neglected; 

Bottcher(SO) # 

Ellipsoids 
multi -phase 

£*  -  c-m 

Short  range  statistical 
variations  neglected;^ 
Polder  and  Van  Santen^* 

1 - 

Ellipsoids 

2-phase 

£-n\~  £},  =■  (£<■  "^h)  $£*  >£*  contains  interaction 

£h  oy.  £<  £*  <  £m  effects;  DeLoor(14). 

. 

s  U  -1 — _j - 

‘ 


7 


Table  2.1  -  continued 
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Neglecting  all  particle  interactions ,  Maxwell^S)  treated  the  problem  of 

anisotropic  ellipsoidal  inclusions  in  an  isotropic  medium  for  magnetic 

(9) 

materials.  For  dielectric  mixtures  similar  equations  are  obtained  . 

,  (14,  30)* 

Rayleigh  showed  that  the  Maxwell  relation  is  a  good  first  order 

approximation  for  cubically  arranged  spherical  and  cylindrical  inclusions 

(14  30) 

in  a  host  material.  A  mixture  equation  developed  by  Wiener  9 

involves  a  form  factor  which,  however,  depends  on  the  dielectric  constants 

of  the  host  material  and  the  mixture  in  a  complex  way.  Wiener  did  not 

(15) 

account  for  particle  interaction  effects  directly.  Lewin  ,  by  applying 

scattering  theory  to  a  material  loaded  with  regularly  spaced  spherical 

particles,  allows  for  particle  interaction  by  defining  an  effective 

dielectric  constant  as  a  complex  function  of  the  particle  size  and  the 

dielectric  constant  of  the  included  particles.  He  specifies  the  particle 

size  to  be  much  smaller  than  the  propagating  wavelength.  Levan’s  formula 

reduces  to  Rayleigh’s  formula  for  ordered  spheres  when  the  effective 

dielectric  constant  of  the  inclusions  is  assumed  equal  to  the  dielectric 

(14  30) 

constant  of  the  material  in  the  inclusions.  Bruggeman  9  extended 

Wiener’s  theory  and  showed  Wiener’s  ’’formzahi”  (form  number)  to  be  a 

function  of  the  dielectric  constant  of  the  host  material  and  of  the  mixture. 

Bruggeman  allowed  for  long-range  particle  interactions  by  the  use  of 

a  limit  procedure.  Niesel^)  arrived  at  Bruggeman’ s  results  by  using 

( 63) 

an  integration  process.  Pearce^  '  modified  Bruggeman’ s  equations  by 
introducing  an  empirical  factor,  k,  which 

*  Both  De  Loor^14)  and  Bottcher^50^  have  given  critical  reviews  of  the 
literature  in  mixed  dielectrics. 
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is  only  dependent  on  the  type  of  mixture  used.  Another  empirical  formula 

(14  39) 

was  given  by  Lichtenecker  *  which  has  no  obvious  correlation  to  the 
other  theoretical  formula  discussed  above, 

Onsager  was  the  first  to  propose  that  only  a  part  of  the  internal 

field  in  a  dielectric  mixture  is  effective  in  directing  a  dipole.  In 

effect,  Onsager  modified  the  average  value  of  the  internal  field  by  a 

reaction  field  which  accounts  for  long-range  interaction  effects  and  hence 

obtained  a  more  accurate  formula  for  a  dipolar  dielectric  mixture. 

Bo’ttcher^^  used  Onsager' s  reaction  field  theory  to  derive  a  mixture 

formula  for  spherical  particles  in  a  continuous  dielectric.  Polder  and  Van 

Santen^^  generalized  this  approach  for  ellipsoidal  shaped  inclusions  of 

more  than  one  material  and  showed  that  Onsager' s  assumptions  about  the 

reaction  field  were  equivalent  to  the  hypothesis  that  "the  mean  value  of 

the  field  in  the  interior  of  the  particle,  of  dielectric  constant  e . 

surrounded  by  a  medium  of  dielectric  constant  £  in  which  a  homogeneous 
J  m 

field  E  is  maintained  at  a  large  distance  from  the  particle,  is  a  good 
av 

approximation  for  the  actual  mean  field  in  this  special  particle". 

A  further  modification  to  the  mixture  relations  of  Polder  and  Van 

(14) 

Santen  was  given  by  De  Loor  who  assigned  an  effective  dielectric  con- 
stant  of  e  to  the  immediate  surroundings  of  an  included  particle.  Even 
though  no  quantitative  solution  is  given  for  the  value  of  £  ,  De  Loor 

clearly  shows  that  the  mixture  dielectric  constant  must  always  lie  within 
the  range  •  The  latter  condition  thus  leads  to  limiting  values 

of  the  dielectric  constant  of  a  mixture. 

(14) 

This  limit  condition  was  verified  by  De  Loor  by  plotting  the  results  2 ^ 
all  the  other  mixture  relations  for  different  materials.  Hasted  and  Shah^ 
applied  this  criterion  to  various  types  of  brick  and  obtained  favourable 
results.  The  lower  limit  is  always  given  by  the  lowest  dielectric  constant. 


■ 

•  •  • 
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Generalized  mixture  equations  can  be  derived  by  the  use  of  vector 
calculus  and  Maxwell’s  equations.  Taylor^*  ^  developed  the  generalized 
equations  for  lossy  anisotropic  dielectric  mixtures  using  the  "average 

t 

field  approximation"  of  Polder  and  Van  Santen.  His  development  which  is 
very  concise  and  formal  is  based  on  the  work  of  Korneenko^^  who  uses 
Gauss’  theorem  to  develop  an  averaging  formula  for  the  parameters  in 
inhomogeneous  media.  Taylor  further  uses  the  contribution  by  Frank^^' 
who  gives  the  penetration  of  a  static  homogeneous  field  in  an  anisotropic 
host  material  into  an  ellipsoidal  inclusion  consisting  of  another  aniso¬ 
tropic  medium.  On  the  development  of  adequate  theories  relating  particle 
interaction  phenomena  in  closed  form,  use  of  Taylor’s  mixture  formulas 
modified  to  allow  short-range  interaction  effects  should  give  a  realistic 
approach  and  solution  to  the  general  problem  of  dielectric  mixtures.  Both 
Kirkwood  and  Brown  developed  statistical  methods  to  solve  the 
particle  interaction  problems.  Brown  shows  that  Maxwell’s  formula  is  a 
second  order  approximation  in  which  statistical  properties  of  particle 
geometry  have  been  neglected.  Bruggeman  and  Bottcher’s  formula  try  to 
account  for  interaction  phenomena  and  agreement  to  the  third  order  with 
Brown's  treatment  requires  assigning  values  to  one  of  Brown’s  statistical 
parameters  which  for  Bruggeman ’s  formula  is  twice  as  large  as  for  Bottcher's 
formula.  To  decide  which  of  these  two  is  the  best  approximation  a  detailed 
calculation  of  the  statistical  particle  distributions  must  be  made. 

Taylor’s  approach  uses  basic  assumptions  about  the  mixture  identical  with 
those  of  Bottcher  and  Polder  and  Van  Santen  and  develops  the  general 
mixture  relations  for  dyadic  permittivities  for  inclusions  with  an  ellip¬ 


soidal  surface. 


} 
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This  chapter  will  develop  Taylor’s  equations v  '  in  modified  form, 

using  De  Loor’s  designation  of  an  effective  permeability  e*  of  the  particle 
(14) 

surroundings  .  The  electrostatic  problem  of  finding  the  ratio  of 
internal  field  to  the  average  field  at  a  large  distance  from  the  particle 
will  be  solved  for  various  practical  particle  geometries.  The  resultant 
equations  will  be  tested  against  experimental  data  for  the  cellulose-air- 
water  mixture  as  found  in  wood  species. 

DERIVATION  OF  AVERAGING  LAW 


It  is  most  important  in  mixture  theory  to  find  the  relationship 
between  the  electric  field  internal  to  an  inclusion  and  the  electric  field 
in  the  surrounding  medium.  To  do  this,  use  is  made  of  Gauss’  Divergence 
Theorem  which  relates  volume  characteristics  of  vector  quantities  to  their 
values  at  the  surface  bounding  the  volume.  Application  of  the  Divergence 
Theorem  to  an  individual  inclusion  in  a  host  medium  and  satisfying  the 
boundary  conditions  at  the  surfaces  between  the  two,  yields  equations 
relating  the  mean  value  of  the  dielectric  constant  in  the  inclusion  to  the 
average  value  of  the  dielectric  constant  in  the  surrounding  medium. 

A  typical  mixture  is  schematically  represented  in  figure  2.1. 


Fig.  2.1  -  Multiphase  Mixture. 


. 

' 


. 
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An  arbitrary  volume  v  is  selected  whose  surface  lies  inside  the 

h 

medium.  Within  this  volume  y  of  surface  y  consider  a  statistically 
large  number  of  inclusions  with  the  volumes  V_^  and  surfaces  ,  (i=l,  2,  .. 
where  the  subscript  i ,  refers  to  a  particular  inclusion.  Let  the  total 
volume  and  surface  be  denoted  respectively^  by 


N 

v=v +X 

h  i=l 


V. 

1 


=  V  +v 

h  I 


and 

N 

S,  =S  +  .11  S. 
h  o  i=l  l 


S  +S_ 
o  I 


(2.1) 


(2.2) 


For  any  vector  X  which  has  continuous  first  partial  derivatives  in  V 
and  on  its  surface  S,  the  divergence  theorem  states 


/yV*  AdV  =  £  n*  AdS 


(2.3) 


where n  is  a  unit  vector  normal  to  the  surface  S  drawn  outward  from  V. 


Consider  now  two  arbitrary  scalars  <J>  and  ip  . 


Setting 

t  =  (pVijj  (2.4) 

and  taking  the  divergence  of  equation  2.4  gives 

V-  (cf>Vij,)  =  V<j>*  V\[j  +  <pV2ip  ^  5) 

Application  of  the  divergence  theorem,  equation  (2.3),  over  the  volume 
and  use  of  equation  (2.5)  yields  the  result 


V  .(Wt)dv-  <t^ds  - 

vh  j\  °n 


(2.6) 


The  negative  sign  appears  because  the  normals  to  the  surfaces  S_^  are 


« 
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directed  into  . 


The  following  quantities  are  now  defined, 

Vip  =  -e>g 

where  e  is  a  dyadic  permittivity  and 


(2.7) 


V<f>  =  u 


(2.8) 


a  unit  vector. 

In  a  source  free  homogeneous  anisotropic  medium  Maxwell’s  equations  yield 
V*  £-2=  0  (2.9) 

Equation  (2.6)  is  now  rewritten  by  the  use  of  equations  (2.7),  (2.8),  and 


(2.9)  as 


N 


'vhH\dv  ■  Jg*e h\ds "  ihjs.*h\ds 

o 


(2.10) 


At  the  boundaries  between  the  host  material  and  the  inclusions  the  boundary 
conditions  are 


e.*E.  =  e*E 

1  xn  h  hn 


(2.11) 


These  boundary  conditions  hold  for  the  general  case  of  a  source  free  material 
with  a  complex  dyadic  permittivity.  Applying  equation  (2.11)  to  the  last 


term  in  equation  (2.10)  one  obtains 


N  I  _  N 


i-l  I  s  *e'hV 

i 


/ 


ds  i=l  J  S. 


6  e  .•  E  .ds 
1  ^n 


(2.12) 


Reapplying  the  divergence  theorem  over  the  volume  V^in  the  same  manner  as 
before,  equation  (2.12)  becomes 


N  r 

L  <pc.'E.  ds 


N  f  +  _  + 

.E  /  u*  £  . •  E . dv 

1=1JV.  1  1 


(2.13) 


' 
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Furthermore,  the  basic  assumption  is  made  that  the  medium  as  a  whole  can  be 
characterized  by  an  equivalent,  average,  anisotropic  dielectric  constant, 

—  — 3^ 

e,  with  a  corresponding  average  uniform  field, E  maintained  at  a  large 

ciV 

distance  from  the  inclusion.  This  is  the  Polder  and  Van  Santen  hypothesis 

which  accounts  for  long-range  interactions,  but  neglects  short-range 

statistical  variations  in  the  particle  geometry.  For  low  value  concentra- 

N 

tions  of  inclusions  ( Jj^A-<l)this  is  indeed  a  good  approximation. 

Hence,  by  definition,  the  following  equation  holds 


-> 

E 


av 


v(Jv  Adv  +  ii  Jv  ,Vv) 

h  ^  l 


(2.14) 


The  first  term  of  the  R  H  S  of  equation  (2.10)  can  now,  by  the  use  of  the 
boundary  conditions,  be  written  in  a  manner  similar  to  equation  (2.13)  as 


-> 


s  *H*h. ds  = 


"n 


„<{>£•  E  ds 
S T  av 
o  n 


=/, 


u*£*Edv 


Vu-^E 


av 


(2.15) 


Substituting  for  the  R  H  S  of  equation  (2.10)  from  equation  (2.15)  and 
(2.13)  respectively  and  noting  the  direction  of u  is  arbitrary,  yields 


V£-E 


av 


.  -*■  N 

e  *E  dv  +  .1 
V,  n  h  i=l 
h 


e  *E  dv 
vi 1 1 


(2.16) 


Exact  calculations  of  the  fields  in  equations  (2.14)  and  (2.16)  is 
very  difficult  in  general  because  of  the  many  complicated  interaction 
phenomena  occur ing  within  any  one  medium.  However,  the  analysis  can  be 
carried  further  if  it  is  assumed  that  each  inclusion  can  be  treated  as  an 


individual  particle  imbedded  in  an  anisotropic  but  homogeneous  material  of 

->  * 

permittivity  e  and  a  field  E^  .  It  is  further  assumed  that  all  inclusions 
are  of  a  shape  which  allows  a  homogeneous  field  to  exist  within  the 


* 


M 


This  is  Polder  and  Van  Santen’ s  hypothesis  and  is  valid  at  least  for  ^ V; « 1 


I 


' 


■ 
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inclusion  when  the  field  in  the  surrounding  host  is  also  homogeneous.  The 
latter  statement  allows  use  of  all  ellipsoidal  inclusions  and  thin  con- 
focal  ellipsoidal  shells 

It  follows,  then,  that  the  field  -within  each  inclusion  is  constant''” 
under  these  assumptions^*  ^*  ^*  61)^  stated  generally 


E.  =  T.E  (2.17) 

l  i  av 


The  quantity  ft  is  a  tensor  of  rank  2  which  in  principle  can  be  determined 
for  each  particle  geometry  satisfying  the  previous  conditions. 


Let  the  average  values  of  K 

h 


and  E . 

l 


be  defined  by 


<E  >  =  ~  I  E  dv 
h  V,  JV,  h 
h  h 


(2.18) 


<E  .  >  =  —  E  dv 
i  V.JV.  i 

l  i 


(2.19) 


Thus  equation  (2.14)  can  now  be  written  as 


->■ 


N 


VE  =  V-<E,>  +  .E_  V .  <E  .  > 
av  h  h  i=l  i  l 


(2.20) 


and  similarly,  for  any  homogeneous  medium,  equation  (2.16)  becomes 


+  ^  N 

VeE  =  V  e  »<E  >  +  .L  V.e:<E.>  (2.21) 

av  h  h  h  i=l  li  i 

Elimination  of  <E,  >  between  equations  (2.20)  and  (2.21)  and  use  of  equation 

h 

(2.17)  finally  yields 

(2.22) 


1  N  _ 

e  “  eh  ~  Vi=l  (£i  “ 


e,  )V.T. 
h  i  i 


The  condition  of  a  constant  internal  field  is  a  direct  result  of  the  finite 
solution  to  Laplace’s  equation  in  ellipsoidal  coordinates. 


'  fc 
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Thus  equation  (2,22)  expresses  the  dielectric  constant, e  ,  of  the 

mixture  in  terms  of  the  dielectric  constants  P  and  r  of  the  inclusion  and 

i  h 

host  material  respectively.  Each  inclusion  contributes  to  the  mixture  in 
proportion  to  its  volume  filling  f actor V./V  and  each  contribution  is 
weighted  by  the  ratio  of  the  internal  inclusion  field  to  the  average  field 
in  the  mixture  as  given  byT.  . 

jL 

The  dyadic  T^  can  be  evaluated  by  methods  of  electrostatics  as 
developed  in  the  next  section. 


■ 
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INTERNAL  FIELD  EVALUATION  FOR  DIFFERENT  SHAPED  INCLUSIONS 

Inclusions  in  many  mixtures  are  small  compared  to  a  wavelength  even 
at  microwave  frequencies.  Hence  the  electromagnetic  solution  to  the  inter¬ 
nal  field  is  quasi-s tat ionary  which  expressed  mathematically  is 

w  «  -°  (2.23) 

a 

Where  c  is  the  velocity  of  light  ^  is  the  characteristic  linear  dimension 
o 

of  the  inclusions. 


Internal  fields  in  inclusions  can  be  solved  in  closed  form  if  the  in¬ 
clusions  are  ellipsoidal.  The  electrostatic  problem  of  an  ellipsoid  in  a 

/  2  ^  a  61  ^ 

homogeneous  field  has  been  treated  by  different  authors^  >  »  » 

However,  an  even  more  general  shaped  inclusion  would  be  in  the  form  of  an 
ellipsoidal  shell.  In  practice,  such  a  shape  is  found  for  example  in  the 
structure  of  wood.  It  is  the  solution  to  the  latter  problem  which  will  be 
derived  in  this  section.  Solutions  for  inclusion  shapes  such  as  spheres, 
spheroids,  needles  and  discs  are  all  contained  in  the  solution  of  the  more 
general  shell  problem. 

A  schematic  representation  of  a  homogeneous  ellipsoidal  shell  in  a 
homogeneous  medium  is  given  in  figure  2.2. 


Fig.  2.2  -  Homogeneous  ellipsoidal  shell  in  homogeneous  medium. 


. 
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In  a  source  free  medium  the  relation  between  the  dielectric  displace- 
ment,  D  ,  and  the  electric  field,  E,  is  given  by 

divD  =  div£*E  =  0  (2.24) 

By  introducing  the  electric  field  potential 
E=  -grad  <f> 

Maxwell's  first  equation 
curlE=0 

is  immediately  satisfied  and  hence  equation  (2.24)  becomes 

div  (E'gradtfO  =  0  (2.25) 

I 

For  homogeneous  media  I  is  a  tensor  of  rank  two  whereby  equation  (2*25) 
reduces  to  Laplace's  equation 

V2<1>  =  0  (2.26) 


In  ellipsoidal  co-ordinates  Laplace's 


equation  becomes 


(2) 


v24>=- 


4 


(€-n)  (C-O  (n-c) 


(n  ?)R£3?(R£3£)+(c  ?)RqaT1(Rri'^')+^"ri)RC3C(Rc'af) 


=o 


(2.27) 


where 

=  / (u+a2)  (u+b2)  (u+c2) 


A  uniform  electric  field  E  is  applied  to  the  material  lying  along  one 


* 

of  the  ellipsoidal  axis  . 


Hence  the  potential  due  to  E  is 


cp  =  -Ex  =  -e/ (£+a2)  (rj+a2)  ( 4+a2) / (b2-a2)  (c2-a2)  (2.28) 

o 


Superposition  can  be  used  to  obtain  the  case  where  E  lies  in  an  arbitrary 
direction. 

From  here  on  it  is  assumed  that  all  three  regions  are  isotropic  but  the 
mixture  itself  may  be  anisotropic.  In  principle,  the  following  procedure 
applies  equally  well  for  anisotropic  media  when  the  formalism  of  tensor 
analysis  is  employed. 


' 


* 


. 
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Outside  the  shell  the  potential  is  modified  due  to  the  presence 
of  the  shell  and  the  enclosed  ellipsoid,  whence 


4>1  =  +<j>' 

I 

where  it  is  assumed 

<f>’  =  4»0F(O 


(2.29) 


Therefore 

cf>1  =  <j>o(l+F(S)}  (2.30) 

^ 1  was  chosen  to  make  the  functional  dependence  of  <j>  on  n  and  £  the  same 
as  for  <j>0  in  order  that  the  boundary  conditions  can  be  satisfied  at  the 
outer  shell  surface  C=0  ,  for  arbitrary  n  and  C  . 


Substituting  equation  (2.29)  into  equation  (2.27)  gives  a  second  order 


differential  equation 


d2F 


dF  d 


df2”  +  dE,  d£  l0S  R^^+a2)  0 


(2.31) 


Two  solutions,  easily  verified  by  substitution,  can  be  found  for  equation 
(2.31),  namely^2,  3) 

F  =  constant  (2.32) 


and 


(2.33) 


F  =  A 


The  upper  limit  in  equation  (2.33)  is  determined  by  the  condition  that  the 
perturbation  effect  due  to  the  ellipsoidal  shell  and  inclusion  must  go  to 


zero  at  infinity. 


In  region  2  the  potential  must  have  the  same  dependence  onq  and  £  as  in 


' 


■ 


region  1,  hence 
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<J>  =  C(p  +  d>'  T 

2  0  Y 

where 

=  *0F'(O 

or 

<P2  =  4>q{C  +  F'  (5)} 

with  F (£),  a  solution  to  equation  (2.32),  given  by 
F'(5)  =B( 


(2.34) 


(2.35) 


In  region  3,  however,  the  behavior  of  the  field  due  to  the  potential 
form  given  by  equation  (2.30)  and  (2.34)  becomes  irregular  at  the  center 
of  the  ellipsoid, /is  £->  - c 2  the  integral  function  F  or  F'  behaves  as  /^+c2 
but  the  potential  gradient  behaves  as  l//^+cz  and  hence  at^=_c2  the  field 
becomes  infinite. 


The  only  remaining  solution  which  satisfies  equation  (2.32)  and  the 
boundary  conditions  on  the  inner  surface  of  the  shell,  for  arbitrary  n 
and  £  is  given  by  equation  (2.32).  Therefore  the  potential  in  region  3  is 
written  as 

d>  =  Dd>  (2.36) 

y3  y0 

A  comparison  of  equations  (2.28)  and  (2.36)  indicates  that  for  a 
uniform  applied  field  t)  the  field  interior  to  the  inner  ellipsoid  is  also 
uniform.  However,  the  field  within  the  shell  is  nonuniform  due  to  the 
added  dependence  on  £,  see  equation  (2.34). 


Continuity  of  the  potential  across  boundary  1-2  requires 


' 

f 
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(2.37) 


Hence,  by  equations  (2.30)  and  (2.34)  and  a  comparison  of  (2.33)  and  (2.35) 
one  obtains 

1  +  F (0)  =  C  +  F'  (0)  (2.38) 


Similarly  at  boundary  2-3  it  is  necessary  that 


S=-k2 


(2.39) 


The  surface  £=-k2  is  the  inner  surface  of  the  shell.  Just  as  the  axes 
of  the  outer  ellipsoid  are  given  by  a,  b,  and  c,  so  the  axes  of  the  inner 
ellipsoid  are  designated  as  a',  b',  and  c'  which  are  defined  as 

a’ 2  =  a2-k2 

b'2  =  b2-k2  (2.40) 

c ' 2  =  c2-k2 


Thus  equation  (2.39)  yields 
C  +  F’ (-k2)  =  D 


(2.41) 


Furthermore,  the  normal  derivatives  must  be  continuous,  hence,  at 
boundary  1-2 


e  i  (' 


i  ■ 


ivhi3£ 


5-0 


=  E 


5=0 


(2.42) 


* 

where 


* 


In  generalized  coordinates  the  element  of  length  is  given  by  du=h2dx. 


(2.43) 


hj  =  Ac-n)(C-C)/2R 


Substitution  for  <j>1  and  §  results  in 


h  fl+^<0)  ]!p 5=o+  <♦„!£> 


5=0 


=  e  [C+  F'(0)]||2- 
2 


4-  (  a  8F\ 

5=0+  W0  3f > 


(2.44) 


?=o 


which  by  the  use  of  equations  (2.28)  and  (2.33)  reduces  to 


* 


e1[l+F(0)]-  =  e2[C+  F’(0)]-2B 


(2.45) 


abc 


At  boundary  2-3 


e  (—  1A 
2  n i 


) 


S=-k2 


e  (-Alii) 


C=-k2 


(2.46) 


which  similarly  reduces  to 

e2[c+  F'<-k2>  ■  sD 


(2.47) 


Equations  (2.38),  (2.41),  (2.45),  and  (2.47)  can  now  be  solved  for 
the  constants  A,  B,  C  and  D  with  the  aid  of  equations (2 . 33) ,  and  (2.35), 

The  ratio  of  the  electric  field  internal  to  the  region  3  to  that  of 
the  field  in  region  1  can  be  expressed  as 

E _  D(fro  =  D 

E1  <j)i  "  1+F(0  (2.48) 


*  3F 

The  term,  — ,  gives  rise  to  the  negative  signs  in  equation  (2.45). 


' 
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* 

At  distances  large  compared  to  the  inclusion  size ;  f(£)-*0  whereby 


la 

El 


E- 


av 


=  D 


(2.49) 


Depolarization  coefficients  for  an  ellipsoid  are  defined  as 


abc  f  d? 

~2  Jq  (C+a^)R5 


(2.50) 


V  Z  rs 

with  n  and  n"  defined  similarly  by  replacing  5  and  a2byrj",  b2  and  r  c2 

respectively.  From  equation  (2.33)  and  (2.50) 

F(0)  =  2nXA  (2.51) 

abc 


For  the  inner  ellipsoid  the  depolarization  coefficients  can  be  redefined  as 


OO 


fx  abc  j  d£ 

2  J_k2.(£+ a2)^ 

(2.52) 

hence 

*<-k2>  “  a'b'c- 

(2.53) 

Use  of  equations  (2.53)  and  (2.51)  in  (2.38), 

(2.41),  (2.45), 

(2.47)  yields 

?  x 

e  =  1  +  ~(A-B) 

(2.54) 

(2.55) 

By  equation  (2.48)  it  is  possible  to  compute  a  correction  factor  to  the 
internal  to  external  field  ratio.  For  low  volume  filling  factors, 
however,  this  correction  is  negligible. 


\ 


'  *  i 


. 
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(2.56) 


e  C  =  e  D  -2B 
2  3 


(n’X-l) 

V'  e2 


(2.57) 


where  V=abc  and  V^a’b’c’ 


Eliminating  C  between  equations  (2.55)  and  (2.57)  yields 


(2.58) 


From  equations  (2.54).  and  (2.56) 


2A{n  (  Ej-  e^)  +  =  2B-^  +  e1  _  e2 


(2.59) 


Substituting  for  B  from  equation  (2.58)  into  (2.55)  and  combining  the 
result  with  equation  (2.54)  gives 


(2.60) 


Substituting  for  A  from  equation  (2.59),  for  B  from  equation  (2.58)  and 
rearranging  the  result  finally  gives 


(2.61) 


Therefore,  the  solution  to  equation  (2.49)  is  given  by  the  value  of  D 
in  equation  (2.61)  which  thus  yields  the  ratio  of  the  field  in  the  inner 
ellipsoid  to  the  average  field  in  the  surrounding  material. 

Similarly  by  equations  (2.30)  and  (2.34)  the  ratio  of  the  field  in  the 
shell  to  the  average  field  in  the  material  is  given  by 


M 
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(E*  /  E^v)  =  C  +  F'(5)  ;  -k2<C<0  (2.62) 

T 

Hence,  the  electric  field  in  the  shell  is  not  uniform  because  of  the 

dependence  on  £ .  For  thin  shells  one  can  average  over  the  thickness  of  the 

shell  to  obtain  a  mean  value  of  the  electric  field  in  the  shell.  To  do 

this  the  value  of  the  integral  function  F(^)  at  the  shell  boundaries  can  be 

found  in  terms  of  the  depolarization  coef f icientsc  Thus  the  sum 
of  F’(f  )  at  the  two  shell  boundaries  is  given  by 


F'U) 


=B 


d? 


-k" 


(£+az)Rc 


dC 


(?+az)R, 


2B,  ,x  V  x. 
— (n  r+n) 


(2.63) 


If  a  linear  variation  of  across  the  shell  is  assumed,  the  average  field 
in  medium  2  is  obtained  by  using  half  the  value  of  equation  (2.63)  in 
equation  (2.62)  which  yields 


av 

Solving  for  C  and  B  from  equations  (2.55)  and  (2.58)  and  using  (2.61)  to 
eliminate  D  yields 


E? 

w 

av 


ei[e2+(£aF2,)(n'*"nX  F)]  ’ 


{[£2+(e2-£3)(nX  rix)][£1+nX(e2-ei)]-nX  ^ 


(2.65) 


Hence,  equations  (2.61)  and  (2.65)  give  the  relationship  between  the 
fields  in  regions  3  and  1  and  regions  2  and  1  respectively.  The  latter 
equation  is  a  good  approximation  for  thin  shells  only. 


Setting  e  =e  in  equation  (2.61)  gives  the  well  known  solution  for  the 
3  2 

internal  field  in  ellipsoidal  inclusions  of  dielectric  constant  s2  dispersed 
in  a  homogeneous  isotropic  medium  of  dielectric  constant  £^ 

-1 


^int 


Ex 

av 


=  e1[£1  +  nX(e2  -  e^l 


(2.66) 


' 

' 

, 


26 


Similarly,  for  a  spherical  shell  of  internal  and  external  radii  b  and  a  and 

dielectric  constant  e3  and  respectively  in  a  host  material  of  dielectric 

* 

constant  e  ,  equation  (2.61)  becomes 


?  =9e2e,[(£2+2el)(2e:2+e3)  •  2&'<e  -e,)e  -e  ) ]  1  (2-62) 

iiaV  z  1  z  a  2321 


,-1 


E 


Equivalent  relations  can  be  written  for  the  y  and  z  directions  of  the 
electric  field  by  replacing  superscripts  x  by  y  or  z. 


However,  the  electric  field  in  the  vicinity  of  an  inclusion  is  not 
exactly  equal  to  that  of  the  average  field  in  the  mixture.  This  is  due  to 
the  spatial  variation  of  the  field  in  the  neighborhood  of  a  particle,  see 
for  example,  equation  (2.48).  In  order  to  account  for  this  variation  one 

k 

can  introduce  an  artificial  dielectric  constant,  8  as  suggested  by 

(14) 

DeLoor  ,  to  replace  the  dielectric  constant  of  the  material  surrounding 
the  inclusions.  Alternatively,  the  correct  field  variation  can  be  calculated 
by  evaluating  the  expression  1+F(£). 

Unfortunately,  both  procedures  require  a  detailed  knowledge  about  the 

statistical  variation  in  particle  location.  In  the  case  of  solids,  relatively 

little  is  known  about  their  statistical  properties  and  hence  the  solution 

to  the  problem  of  finding  the  exact  depolarizing  fields  in  dielectric 

mixtures  can  only  be  estimated  by  empirical  methods.  Nevertheless  bounding 

k 

values  can  be  determined  for  the  effective  dielectric  constant  e  which,  as 


* 


Landau  and  Lifshitz  solved  this  problem  (see  reference  2,  p.45)  and 
gave  the  solution  for  the  sperial  case  s1=s2=l  (vacuo)  and  e=  e. 


. 


. 


/ 


27 


De  Loor  showed  for  non-dissipative  media,  always  falls  between  the  value  of 
the  host  and  mixture  dielectric  constants. 


In  actual  mixtures  one  often  finds  a  thin  shell  of  included  material, 

such  as  absorption  layers,  inside  a  host  medium  in  which  case  e  =e  . 

1  3 

k 

Maying  the  substitution  e  ^—63=  e  in  equation  (2,65)  and  substituting  the 
results  for  1\  in  equation  (2.22)  finally  yields 


=  v  (e  )e  [e  + 


*  TT  T 

e  -e?  ,  ,  V\  , 

- ~(n  !-n  — )  ] 


'h  VV  V"  L °2  2 

VT 


j  j  V' 


{[e2+(e2-£X)(nj~  -  nj  )  ]  [e*-hu  (e^e*)  (e2-e*j  } 


. 


(2.68) 


where  is  the  normalized  filling  factor  and  the  subscript  j  refers  to 
the  coordinate  axis. 

Equation  (2.68)  relates  the  mixture  dielectric  constant, to  the 
dielectric  constant  of  the  host  material, ,  with  thin  ellipsoidal 
inclusions  of  dielectric  constant  c2  anc*  volume  filling  factor  v2« 

Limiting  values  for  the  mixture  dielectric  constant  are  obtained  by  solving 

3*?  * 

equation  (2,68)  with  e  ®  ehand  e  If  spherical  and  ellipsoidal  inclusions 

also  exist  within  the  mixture  then  other  terms  will  appear  in  the  equation 
as  a  direct  result  of  the  substitution  for  the  various  values  of  T.  in 

l 

equation  (2.22). 


Since  the  boundary  conditions  to  the  above  problem  were  stated  without 
any  restriction  on  the  dielectric  constants,  the  solutions  are  valid  for  lossless 
or  lossy  materials.  The  only  restriction  is  that  the  assumption  about  the 
uniformity  of  fields  be  applicable.  In  order  to  calculate  the  low  and  high 
frequency  dielectric  constants  of  a  mixture,  the  appropriate  dielectric 
constants  for  the  included  and  host  materials  are  used  in  the  equations. 
Generally,  the  equations  developed  here  will  be  valid  for  any  frequency  and 


I 
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temperature  if  the  corresponding  values  of  the  individual  dielectric 
constants  are  used. 

Practical  Significance  of  Thin  Shell  Inclusions 

It  can  readily  be  shown  that  for  the  case  of  a  spherical  thin  shell 
with  b-a,  e3=ei’  anc*  £2>>el  t^ie  field  internal  to  the  shell  is  approx¬ 

imately  2/3  times  the  average  field  in  the  mixture.  Ellipsoidal  shells 
would  tend  to  enhance  this  internal  field  even  further  in  the  direction  of 
the  major  axis.  This  internal  field  for  a  shell  inclusion  is  about  twice 
as  large  as  the  internal  field  found  for  spherical  inclusions  under  the 
same  assumptions.  Increasing  the  internal  field  causes  an  increase  in  the 
dielectric  constant  of  the  mixture.  This  is  significant  since  present 
mixture  theories  applied  to  materials  assumed  to  have  spherical  or 
ellipsoidal  inclusions  yield  values  for  the  dielectric  constant  which  are 

k 

lower  than  the  experimental  values. 

It  is  suggested  by  the  author  that  the  neglect  of  thin  shell  inclusions, 
when  they  exist,  will  be  an  important  contributing  factor  to  such  discrepan¬ 
cies.  In  any  case,  it  becomes  apparent  that  the  geometry  of  the  inclusions 
is  an  extremely  important  factor  in  dielectric  mixture  calculations. 


VALIDITY  OF  THE  MIXTURE  RELATION 

In  the  derivation  of  the  general  mixture  relation  the  following 
assumptions  were  made:  1.  The  mixture  can  be  described  by  an  average 
dielectric  constant,  c,  and  an  average  electric  field,  Eav .  2.  Inclusions 


(62) 


Hasted  and  Shahvvw  found  that  the  limiting  values  of  mixture  dielectric 
constants  given  by  De  Loor  (14)  fall  below  the  experimental  values  found 
for  brick.  It  is  very  likely  that,  in  brick  thin  shells  of  water  are 
formed  along  the  inside  walls  of  air  cavities  and  around  particles  making 
it  necessary  to  use  the  internal  field  predicted  by  equation  (2.65) 
rather  than  equation  (2.66)  commonly  used. 


. 
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are  treated  individually  as  being  present  in  a  material  with  the  properties 
just  stated  under  (1) .  Hence,  particle  interaction  effects  are  neglected 
insofar  as  the  value  of  the  calculated  internal  field  accounts  only  for 

•k  _  ' 

the  individual  reaction  fields  of  each  inclusion  and  e is  an  average 
dielectric  constant  of  the  remaining  mixture.  3.  Quasistationary 
arguments  must  be  applicable  to  the  mixed  system. 


Assumption  (1)  is  a  classification  of  the  problem  and  its  solution 
rather  than  a  restriction  on  the  theory.  It  does,  however,  require  the 
validity  of  assumption  (3)  because,  for  inclusions  comparable  in  size  to 
a  wavelength,  the  field  cannot  be  uniform  any  more  and  hence  the  mixture 
cannot  fairly  be  described  by  an  average  field  distribution.  It  has  been 
shown  that  inhomogeneity  of  the  local  field  in  the  neighborhood  of  a 
small  particle  makes  no  primary  contribution  to  the  average  field  (see 
Kirkwood  ^^  )  .  Hence  small  irregularities  in  the  shape  of  the  particles 
can  be  neglected. 


Validity  of  assumption  (2)  requires  that  the  square  of  the  normalized 

N 

volume  filling  factor V^/V,  is  small  compared  to  unity.  For  large 
filling  factors,  deviation  from  the  theoretical  behavior  is  to  be 
expected  even  though  the  interaction  phenomena  is  often  a  second  or 


third  order  effect. 


(21) 


This  author  is  of  the  opinion  that  a  large  number 


When  a  particle  is  placed  inside  a  cavity,  the  same  size  as  the 
particle,  in  a  dielectric  medium,  there  is  a  contribution  to  the 
total  average  field  in  the  medium  due  to  the  polarization  of  the 
particle.  This  added  term  to  the  total  field  is  often  called  the 
’reaction  field’  (10,  11,  12,  14).  In  the  previous  section  the 
reaction  field  was  given  by  the  potential  function  cJm. 

Ey  virtue  of  equation  (2.30)  the  electric  field  next  to  an  inclusion 
is  not  uniform  anymore  but  its  spatial  dependence  is  calculable. 


.  ...  > 

4 

' 


••  ‘ 


, 


. 


50 

of  the  errors  associated  with  the  various  mixture  theori.es  are  due  to  the 


fact  that  mixtures  to  which  the  relations  are  applied  do  not  meet  the  basic 
assumptions  of  the  theory  aside  from  particle  interaction  phenomena.  A 
distribution  of  a  large  number  of  particle  shapes  in  a  mixture  can  contri¬ 
bute  major  errors  if  the  theoretical  relation  simply  assumes  particles  of 
one  specific  geometry.  The  present  generalized  treatment  of  the  problem 
should  yield  more  accurate  estimates  of  practical  mixtures,  eventhough  the 
computational  difficulty  is  increased.  The  latter  objection  is  largely 

negated  by  the  increased  availability  of  direct  access  computer  terminals. 

(2 A) 

All  the  calculations  in  this  thesis  were  done  by  using  APL  programs 
on  an  IBM  360/67  computer.  The  theory  of  this  chapter  is  applied  tc  the 
cellulose- air ~Y/ater  system  in  chapter  6  and  the  results  compared  to  ex¬ 
perimental  data  on  such  a  mixture. 


CHAPTER  3 
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TEMPERATURE  AND  FREQUENCY  DEPENDENCE  OF  DIELECTRICS 

INTRODUCTION 

Knowledge  about  the  temperature  and  frequency  behavior  of  dielectric 
materials  is  essential  for  designing  systems  in  which  dielectrics  play  an 
important  role.  For  example,  in  microwave  power  applicators  the  dielectric 
properties  of  the  treated  material  determine  the  heating  rate,  drying  rate, 
stability,  maximum  attainable  temperature,  energy  distribution  and  electro¬ 
magnetic  coupling  efficiency. 

This  chapter,  divided  into  two  sections,  will  give  the  fundamental 
equations  relating  the  dynamic  properties  of  dielectrics  to  some  basic 
material  constants.  Experimental  verification  of  a  few  of  the  results  is 
given  here  with  the  more  extensive  experimentally  obtained  data  analyzed 
in  chapter  5. 

In  the  first  section  the  treatment  covers  a  generalized  formulation 
of  the  temperature  dependence  of  wide  dispersion,  low  loss  dielectrics. 

The  second  section  gives  a  detailed  analysis  of  the  temperature  and 
frequency  dependence  for  materials  exhibiting  a  Gaussian  distribution  of 
relaxation  times.  Solutions  of  the  integral  equations  are  obtained  by  the 
use  of  numerical  techniques  on  a  digital  computer. 

A.  TEMPERATURE  BEHAVIOR  OF  WIPE-DISPERSION  DIELECTRICS 

/oo  \ 

The  work  of  Gevers  and  Du  Pre  v  J  (1945)  represents  a  major  contri¬ 
bution  to  the  understanding  of  the  temperature  dependence  of  materials 
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which  they  classified  as  "technical  solids"  (low  loss  tangents).  Their 
analysis  shows  that  in  the  case  of  dielectrics  having  either  a  small  thermal 
expansion  coefficient,  or  a  low  value  of  the  dielectric  constant,  the  ratio, 
A,  between  the  temperature  coefficient  of  the  dielectric,  constant  and  the 
loss  tangent  is  approximately  independent  of  the  nature  of  the  dielectric  in 
the  frequency  range  of  1000  Hz  to  10  MHz  and  at  temperatures  between  -50  and 
150  C.  They  show  A  to  have  a  value  lying  between  .04  and  .06. 

Underlying  their  analysis  is  the  basic  assumption  that  relaxation 

phenomena  in  dielectrics  can  be  described  in  a  manner  similar  to  a  chemical 

rate  process  and  that  a  distribution  of  relaxation  times  exists  in  the 

(39) 

material.  Gevers’  work  includes  a  comprehensive  review  of  the  litera- 

(41  44  66  ) 

ture.  His  results  have  been  used  extensively  by  other  authors.  ’  ’ 

(41) 

Garton  has  suggested  that  a  more  generalized  approach  should  be 

taken  to  the  derivation  of  Gevers’  and  Du  Pre's  law.  This  general  law, 
relating  the  ratio  of  dielectric  temperature  coefficient  and  loss  tangent  is 
derived  for  the  first  time  in  this  thesis.  Experimental  data  will  show  the 
new  law  to  be  more  accurate  and  satisfactory  in  practice. 

DERIVATION  OE  A  GENERAL  LAW  FOR  WIDE  DISPERSION  DIELECTRICS 

Any  dielectric  polarization  mechanism  which  is  linear  with  regard  to 

/^0  ^3  00  \ 

the  applied  field  can  be  represented  by  ’  ’ 


oo 


■f(T)  dr 

I  +U2T2 


(3.1) 


V. 


(3.2) 


. 
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where  C  is  the  relaxation  time  of  a  particular  polarization  mechanism,  CJ 
is  the  radial  frequency,  -f  (T)  is  a  distribution  function  of  relaxation 
times,  ^  is  the  value  of  the  dielectric  constant  at  a  frequency  high 
enough  to  make  the  dispersion  of  £’  due  to  the  assumed  polarization 
mechanism  negligible,  £0  is  the  value  of  £'  at  a  frequency  much  smaller 
than  the  frequency  at  which  maximum  loss  occurs,  and  Tg  is  the  shortest 

k 

time  constant  present  in  the  substance. 


It  is  desirable  to  find  a  distribution  function  T  (T)  such  that  the 
frequency  variation  of  the  loss  tangent  is  nearly  constant  over  a  wide 
range  of  frequencies  since  many  solid  dielectrics  exhibit  such  a  behavior. 
By  inserting  a  few  trial  functions  into  equation  (3.2)  it  is  found  that  a 
form  of  f  (T)  which  leads  to  a  constant  value  of  £"  over  a  wide  frequency 
range  is  given  by 


f(T)=  l/T 


(3.3) 


Substitution  of  equation  (3.3)  into  equation  (3.2)  and  using  the 
change  of  variable  CJF=  tan©  under  the  condition  that 


c  <  cjr0  «  | 


(3.4) 


whence 


tan"1  [  WtJ 

gives 


and 


u 


(3.5) 


£"=tt[e-^)/2  (3.6) 


To  is  the  time  required  for  molecules  to  turn  from  one  equilibrium  position 
to  another  multiplied  by  a  probability  factor  which  is  related  to  the  temp¬ 
erature  dependence  of  the  energy  barrier  separating  two  equilibrium 
positions.  (39^  66)  (see  also  section  B  of  this  chapter.) 


From  equation  (3.6)  it  follows  that 
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£0-£„  -  ?-  £''/tc  —  Z  £' -tan  6 /tt 


which  when  combined  with  equation  (3.5)  yields 


2.  £'  Tdn  £ 


Ln 


U)TC 


(3.7) 


(3.8) 


It  is  usually  assumed  that  the  high-frequency  dielectric  constant,  £»> , 

is  independent  of  temperature.  That  this  assumption  is  in  fact  correct  is 

(18) 

shown  by  Van  Vleck  in  his  quantum  mechanical  treatment  of  the  dielectric 

■k  ,j^ 

constant.  Furthermore,  t0  s  the  shortest  time  constant  present,  is  assumed 
to  be  independent  of  temperature  as  well  and  can  be  regarded  as  a  constant 
of  the  material. 


Differentiating  equation  (3.8)  w.r.t.  the  temperature, T  ,  yields 


dfj 

df 


/0  I  „  .  if  ds‘  i 

1  w  i;i“n  Jr  i 


I  }  r  I  4  p  /“ 

io.no  I  •+-  E  — 


/ 


a 


or 


i  a  £' 

£'ian<f  cTt 


2, 

7T 


Ln  ( i  ^ p)  ciT  ^ ^  1  *■' 

"TTan  £  Ln  (T/ToToT 


o 
c - 

rr 


(3.9) 


(3.10) 


When  either 
constant  of 

df^JclT,  may 


the  linear  coefficient  of  thermal  expansion,  or  the  dielectric 

the  material  is  very  large  the  temperature  expansion  term, 

•ki c-,V 

be  appreciable  but  can  readily  be  calculated  and  included 


*** 


There  is  a  temperature  dependence  of  the  high  frequency  dielectric  constant 
only  insofar  as  the  density  of  the  material  changes  with  varying  temperatures 
Experimental  evidence  to  justify  this  assumption  is  given  later  in  this 

chapter.  /(38) 

Based  on  Claus ius-Mosso t ti ’ s  Law,  Gevers  and  DuPre  derived  the  following 

expression  for  the  temperature  dependence  of  £.cc.  ; 

i  cl  C. 


r- 1 

t 


dT 


—  —  oC, 

i 


4-  8c0  -  2,/ty, 


where  c.^  is  the  linear  expansion  coefficient  of  the  dielectric. 


. 


* 


' 
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in  equation  (3.10).  The  temperature  dependence  of  the  loss  tangent  for 

>fC 

the  class  of  materials  treated  here  can  be  written  as 

tons  ck  Tne'a/T  (3.11) 

(3.12) 

Thus,  equation  (3.10)  can  be  rewritten  as 


Hence 


cL  Lh  tan  <S 
clT 


tan 


0 


d  tan  £ 

cPF 


T2- 


A'=  — _  =  &Ln(\/W. J'£  +  &) 

Z'tanS  dT  I  —  —  Ln  ( I /a)Tc)van  £ 


(3.13) 


The  term  in  the  denominator  on  the  R.H.S.  of  equation  (3.13)  is  only 
close  to  unity  when 

_L  Lnfl/urJ-tan  5  «  I 

tt  1  (3.14) 

In  the  special  case  when  n*l,  a=o  and  the  condition  (3.14)  holds, 
the  general  law  given  by  equation  (3.13)  reduces  to  Gevers'  and  Du  Pre's 
result,  namely. 


A  -  j  d£'=  Ju  Ln(  I'/WToj  (3.15) 

t‘ian  s  dT 

The  relation  in  equation  (3.15)  is  a  good  approximation  only  when 
the  loss  tangent  depends  linearly  on  temperature  and  obeys  condition  (3.14) 


Equation  (3.11)  is  based  on  a  theoretical  treatment  given  by  Gar  ton 
and  relates  favorably  to  experimental  results  on  low  loss  materials. 


: 


ratio  of  temperature  coefficient  to  loss  tangent 

©  •  *  ® 
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A  plot  of  equation  (3.13)  and  its  special  case,  equation  (3.15), 


for  an  arbitrary  value  of  Ln(1  /WT^)  is  given  in  figure  3.1  with  n  chosen 


as  unity  and  a  given  as  a  parameter.  The  plot  shows  very  distinctly  the 


expression. 

For  the  sake  of  simplicity  the  denominator  of  equation  (3.13)  was  set  to 
unity  to  obtain  figure  3.1.  The  added  effect  of  the  denominator  is  an 
increased  magnitude  of  the  ordinate  for  nonzero  a-values  with  the  largest 
positive  correction  occurring  at  the  highest  temperatures. 

Figures  3.2  through  3.4  give  plots  of  A  as  a  function  of  temp¬ 
erature  for  mycalex,  porcelain,  and  ebonite  respectively.  The  data  were 
originally  plotted  by  Gevers  (see  part  V  of  Reference  39)  and  compared 
with  a  plot  of  equation  (3.15)  which  gave  rise  to  large  deviations  at  the 
low  temperature  end.  Here  the  general  equation  (3.13)  and  its  special 
case,  equation  (3.15),  are  plotted  and  compared  with  actual  experimental 
data,  quoted  from  Gevers*  work,  indicating  an  improved  correspond¬ 
ence  between  theory  and  experiment,  especially  at  the  lower  temperatures. 
It  is  further  evident  that  the  generalization  of  the  temperature  behavior 
of  the  loss  tangent  still  allows  the  value  of  A  to  remain  fairly  constant 
for  the  class  of  materials  under  discussion.  Values  of  the  exponent  a 
were  chosen  to  give  a  reasonable  fit.  between  the  theoretical  and  experi¬ 
mental  curves. 

(39) 

In  part  III  of  Gevers'  paper  it  is  stated  that  the  ratio  of  the 
temperature  coefficient  and  loss  tangent  for  a  mixture  having  components 
with  both  positive  and  negative  temperature  coefficients  cannot  be  a  con¬ 


stant  anymore.  However, 


vs.  temperature  for 
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* 
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Fig.  3*2  -  Theoretical  and  experimental  plot  of  ( d£ J/dT)/£ 1  tan  £ 
versus  temperature  for  Mycalex. 


ratio  of  temperature  coefficient  to  loss  tangent 
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Experimental 
Equation  (3*15) 

^-Equation  (3*13)»  a=-50  ,10=16 
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Fig.  3*3 


^  ‘  “  2‘00  '  360  400 

Temperature  °K 

-  Theoretical  and  experimental  plot  of  (d£ */dT)/£ *  tan  <f 
versus  temperature  for  Porcelain. 
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100  . "  " . '  . '  200  p  500 

Temper sirre  °K 

Fig.  3.4  -  Theoretical  and  experimental  plot  of  (d£*/dT)/f*tan S  versus  temperate* 
for  Ebonite. 
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Douglas  Fir  at  a  low  percentage  moisture  content  indicates  values  for  A 
falling  within  the  range  of  ,02  to  .08,  (see  figure  3,5)  over  a  narrow  tempera¬ 
ture  range.  Even  though  cellulose  and  water  have  temperature  coefficients  of 
opposite  sign,  a  cellulose-air-water  mixture  such  as  wood  gives  rise  to  a  wide 
distribution  of  relaxation  times,  flattening  the  loss  vs.  frequency  curve, 
and  hence  approximately  obeying  the  basic  assumption  used  to  derive  equation 
(3,13), 


It  will  be  shown  in  chapter  5  that  the  distribution  of  relaxation  times 
exhibited  by  wood  obeys  a  Gaussian  distribution  and  not  a  l/T  distribution  as 
assumed  here.  For  this  reason  the  temperature  variation  of  A  as  a  function 
of  T  calculated  for  Douglas  Fir  by  means  of  equation  (3.13)  is  not  as  great 
as  the  experimental  temperature  variation  which  is  evident  from  figure  3.5. 

It  is  concluded  that  the  ratio  of  the  temperature  coefficient  of  the 


dielectric  constant  to  its  loss  tangent  is  reasonably  constant  for  a  wide 
range  of  dielectric  materials  which  have  an  extremely  wide  dispersion  range 
The  frequencies  at  which  this  law  holds  are  very  much  lower  than  the 
frequency  corresponding  to  the  shortest  relaxation  time  in  the  material. 


_B .  GENERAL  TEMPERATURE  AND  FREQUENCY  DEPENDENCE  OF  DIELECTRICS 

The  analysis  given  in  this  section  is  based  on.  fundamental  equations 
in  dielectric  theory.  By  assuming  a  particular  distribution  of  relaxation 
times  these  equations  are  developed  to  give  the  frequency  and  temperature 
behavior  of  the  material  in  terms  of  basic  material  constants.  In  chapter 
5  the  theory  is  applied  to  a  very  practical  and  naturally  occurring 


** 

* 


, 


ratio  of  temperature  coefficient  to  loss  tangent 
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Equation  (3.13) 
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Temperature  °K 

Fig.  3.5  -  Theoretical  and  experimental  plot  of  (dE’/dT) £* tan G 
versus  temperature  for  Douglas  Fir  at  1%  moisture 
content. 
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dielectric  mixture  and  the  results  compared  with  experimental  data.  It  is 
assumed  tnroughout  the  analysis  that  the  dielectric  dispersion  is  caused 
by  orientation  polarization.  Since  application  of  the  theory  is  to  be 
made  at  microwave  frequencies,  conduction  effects  are  small  and  therefore 
ignored.  Neglect  of  the  atomic  and  electronic  polarization  mechanisms, 
giving  rise  to  a  resonance  type  of  behavior  at  frequencies  in  the  infrared 
range,  does  not  effect  the  generality  of  the  present  analysis. 

Discussions  on  relaxation  times,  activation  energies,  and  their 
dependence  on  temperature  can  be  found  in  the  literature  (References 
19,  38,  39,  40,  44,  45,  47,  49  and  65  ) .  An  effort  is  made  to  combine 
the  work  of  various  authors,  in  as  far  as  it  aids  the  analysis,  to  extend 
the  results  and  to  apply  them  to  the  cellulose-air-water  system. 

For  dilute  solutions  of  dipoles  in  liquids  and  solids  the  Debye 
equations  predict  the  dielectric  behavior  very  satisfactorily.  However, 
these  equations  give  rise  to  only  a  single  relaxation  time.  In  practice 
it  is  found  that  mixtures  of  materials  contain  more  than  one  relaxation 
time.  In  fact,  very  wide  distributions  of  relaxation  times  may  exist. 

Equations  (3.1)  and  (3.2)  of  the  previous  section  have  been  formulated 
to  allow  for  such  a  distribution.  In  the  following  analysis  it  is  assumed 
that  the  spread  of  relaxation  times  obeys  a  Gaussian  distribution  function. 
Such  a  distribution  function  allows  relaxation  times  with  values  between 
zero  and  infinity  and  is  mathematically  stated  as 


(3.16) 


Experimentally,  it  is  observed  that  the  dielectric  dispersion  and  loss 
is  symmetric  on  a  logarithmic  irequency  scale,  liquation  (5,lo)  tall 
behave  similarly  by  making  the  transiormation 


X=  In  (T/Tj 


(3.17) 


■ 


. 

> 


, 
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in  which  l  is  the  relaxation  time,  is  the  most  probable  relaxation  time, 
and  b  is  a  probability  density  parameter.  Equation  (3.16)  substituted 
into  equation  (3.2)  results  in 


£"  =  b(£0-£J 

Vtt 


x0=U  (TJx) 


i 


-PX*  _  , 

& _ ..Cy>.T_  d  X 

v  ,  .  ,2/7^2. 


J  +  to*T 


(3.18) 


The  temperature  dependence  of  £"  is  implicit  in  the  temperature  behavior 

of  the  relaxation  time T .  Following  closely  the  arguments  set  forth  by 
* 

Gevers  the  polarization  mechanism  is  explained  in  subsequent  paragraphs. 

Particles,  such  as  dipoles  and  ions,  give  rise  to  dielectric  losses 
and  are  subject  to  various  force  fields  (Van  der  Waalfs  forces 5  Coulomb 
forces)  due  to  their  surroundings.  Each  individual  particle  thus  has  a 
certain  potential  energy  w.r.t.  its  neighbors.  Cellulose,  the  dielectric 
under  study,  consists  of  long  molecular  chains  which. when  part  of  a  solid, 
amorphous-crystalline  structure  like  wood,  has  a  great  variety  of  orienta¬ 
tions.  Hence  the  assumption  is  made  that  a  large  spread  of  energy  levels 
exists.  A  schematic  presentation  of  such  an  energy  system  is  given  in 
figure  3.6. 

When  the  material  is  in  a  state  of  equilibrium  the  particles 
oscillate,  due  to  Brownian  motion,  in  the  potential  energy  wells.  From 
figure  3.6  it  follows  that  in  order  for  a  particle  to  jump  to  another 
energy  level,  say  it  must  receive  a  boost  in  energy  equivalent  to  the 

Gevers (part  III).  Garton^^  also  treats  this  subject,  but  restricts 
his  analysis  to  the  specific  case  of  materials  with  a  low  but  constant 
loss  tangent,  whereas  Gevers1  arguments  are  more  widely  applicable  and  are 
substantiated  by  others  (44,  45,  46,  66)  who  treat  relaxation  in  a  manner 
similar  to  a  chemical  rate  process. 
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Deviation  from  Equilibrium 


Fig,  3,6  -  Schematic  presentation  of  an  energy  position 
diagram  for  a  typical  solid  dielectric. 


barrier  height,  ,  which  can  occur  through  Brownian  motion  or  by 
collision  with  other  particles  arriving  or  present  at  that  particular 


is  the  absolute  temperature  and  k  is  Boltzman's  constant.  Thus,  the 
number  of  oscillations  about  the  equilibrium  position  which  a  particle 


during  which  a  particle  is  in  its  equilibrium  position  is 


(3.19) 


where  X0  is  the  proportionality  factor. 


*  T ,  consists  of  the  product  of  two  terms.  One  term  represents  the  time 
of  oscillation  between  equilibrium  positions  and  the  second  term  is  a 
probability  factor  related  to  the  population  density  in  the  higher  and 
lower  energy  states.  T0  is  usually  found  to  be  shorter  than  the  time 
required  for  a  molecule  to  change  equilibrium  positions. 
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Hereafter  it  is  assumed  that  equation  (3.19)  describes  the  temperature 
dependence  of  the  relaxation  time  of  the  molecular  polarization  because 
when  an  electric  field  is  applied  to  the  material  the  dipoles  rotate  and 
the  ions  move  in  the  direction  of  the  applied  field  requiring  for  this 
process  the  activation  energy  q. 


The  proportionality  factor  T o  active  in  the  process  is  assumed  inde¬ 
pendent  of  the  activation  energy  and  the  temperature.  In  general  the 

activation  energy  will  depend  on  temperature  as  is  the  case  for  cellulose 
(47) 


mixtures 


If  such  a  dependence  is  approximated  linearly  as 


V=f-  aT 


(3.20) 


where  a  is  a  constant,  then  equation  (3.19)  can  be  rewritten  as 

r  =  r/e^'aT)AT=  T/e  a/*  e  * /kr 


(3.21) 


hence 


_a/k 


(3.22) 


To  e' 

Using  Q  as  the  activation  energy  in  calories  per  mole  of  material,  the 


expression  (3.21)  becomes 

Q/RT 

T=  T0  e 


(3.23) 


where  R  is  the  gas  constant  and  is  1.987  cal/mole  -  K. 


From  equation  (3.17) 


then 


X  =  Ln  T/T 


t-  t  £* 


(3.17) 


(3.24) 


Equation  (3.24)  can  now  be  substituted  into  equation  (3.18)  to  give 

zM  I 
■  W 


£"=  bfgn-g„)  I  d  x 

—  1  i  +(ojrex)2- 


(3.25) 


*0 


where 
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r=r 


Q/RT 

e 


(3.26) 


Hence  equation  (3.25)  and  (3.26)  describe  the  functional  relationship 
between  dielectric  loss,  temperature,  and  frequency,  under  the  assump¬ 
tion  that  polarization  occurs  because  of  molecular  displacement  and  ro¬ 
tation  subject  to  a  Gaussian  distribution  of  relaxation  times. 


In  a  completely  analogous  manner,  the  temperature  and  frequency 


dependence  of  the  real  part  of  the  dielectric  constant  is  derived.  Use 


of  equations  (3.16)  and  (3.24)  in  equation  (3.1)  yields 

e  _ rlX 


6  —  6^  -f  b  (~o  £<*>  ) 

Vff 


(3.27) 


)  +(core*J2r 

u 

To  find  the  temperature  at  which  the  maximum  loss  occurs,  equation  (3.25) 
is  differentiated  w.r.t.  temperature.  Therefore,  under  the  assumption 


that  the  integral  is  continuous , 

roo  >2  z 

d£"=  b  (£„-£„)  /  e 

dT  —7W  7v 


wtex 


d  •  _ 

dT  i+fwred 


TT2T 


d* 


(3.28) 


Carrying  out  the  differentiation  by  the  use  of 

de"=  b£_“  ar 

dT  i>T  cT 

equation  (3.28)  becomes 

d£"--Qb(£0-£j  /  i-klx2cofexLl-(^tey]rix  (3.29) 

dT“  Hr *W~Ji  [i  +  («reTT  di 

xo 

where  f(x  )  is  the  term  due  to  the  lower  limit.  This  term  can  be  neglec— 
o 

ted  if  _ 

(ATmJ  >  fWd*!. 
dT  dT 

where  AT  is  the  maximum  uncertainty  in  the  temperature  at  which  the 
max 

loss  maximum  occurs  (£&  1°C) .  It  can  be  shown  that  the  right  hand  side 
of  the  inequality  is  extremely  small,  namely,  <1  10  ,  over  the  frequency 


and  temperature  range  of  interest  so  that  the  inequality  nolds# 


•  4  -  * 
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Similarly,  differentiating  equation  (3.27)  w.r.t.  temperature  gives 

-oo 


d£  =  +  <?b(£e-£.)  A  (core)  zch 

dT  dr  RT^I/rr  JXa  [|  +  (wre*)2']1 


(3.30) 


term  ^  0 


In  order  to  find  the  maximum  loss  temperature  it  is  required  that 


dj"~ 

dT 


—  o 


ifT=T  and  U)—  constant 


since  the  point  of  maximum  loss  occurs  at  T=T  by  definition.  The 
right  hand  side  of  equation  (3.29)  vanishes  if  and  only  if 


u)Tey=  Ca/F=  i 


(3.31) 


Hence,  at  T=T  equation  (3.31)  can  be  written  as 
_*  Q/R~Tm 

cor  =  CJTn  e  =  I 


(3.32) 


where  is  the  temperature  at  which  the  maximum  loss  occurs. 

In  order  to  evaluate  gome  of  the  constants  in  equations  (3.25)  and 

(3.27),  condition  (3.32)  is  substituted  in  these  equations  so  that 


£"(TJ=  b(L0-U  /  e-^-e*.  dx 

/  i  +  e 


(3.33) 


Let 


■°°  ,  Z.2. 

/  —  k  X  X  i 

QM  =  /  e  _ dx 


(3.34) 


and  therefore 


// , 


(£.-£-)  A  =  £ 

Vrr  g (x) 


(3.35) 


\ 


/ 
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Also, 


(3.36) 


Now  let 


f  fx)  = 


(3.37) 


Hence 


£„-£xb=  £  '(Tm)  -  <f, 

VW  f  (x) 


(3.38) 


Substitution  of  equation  (3.35)  into  equation  (3.38)  gives  the  value  for 


£»  -  £'(V  -  4'fU  im 


(3.39) 


3W 


Using  equations  (3 . 35) ,  (3.38)  and  (3.39)  in  equations  (3.25), 
(3.27),  (3.29)  and  (3.30)  finally  gives  the  complete  macroscopic  di¬ 
electric  properties  of  a  material  w.r.t.  frequency  and  temperature  in 
terms  of  a  few  basic  material  constants.  These  constants,  the  shortest 


relaxation  timeT^  ,  the  activation  energy  of  polarization,  Q,  and  the 
most  probable  relaxation  time  "l  are  related  through  equation  (3.32). 


The  dispersion  factor  b  in  the  assumed  Gaussian  distribution  is  deter¬ 
mined  by  fitting  the  theoretical  data  to  the  experimental  data. 

EVALUATION  OF  THE  BASIC  CONSTANTS 

Activation  energy  is  the  energy  required  to  make  a  particle  jump 

(47) 


from  one  equilibrium  position  to  another.  The  value  of  this  energy 
barrier  can  be  obtained  from  the  slope  of  a  plot  of  LnTT  versus  l/T 
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ig .  3#7  -  Log(TT)  versus  1/T  for  Pinewood  at  various 
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For  the  case  of  pinewood  the  variation  oft  vs  T  is  given  in  the  litera- 
ture^^  .  The  data  are  plotted  with  moisture  content  as  a  parameter  in 
figure  3.7,  The  activation  energy  can  be  shown  to  obey  the  relation 
<?=  RtLnKTj-Ln^r.^/LO/T,)  "(l/T,.)]  (3.40) 

However,  the  deviation  of  the  points  in  figure  3.7  doesn't  allow 
for  a  very  accurate  determination  of  Q  vs.  moisture  content.  Hence,  an 
iterative  numerical  procedure  was  used  to  evaluate  Q  using  equation  (3.2) 
which  can  be  written  as 

<?=  RT  Ln  (T/T„)  (3-41) 

as  pointed  out  in  the  previous  section,  the  temperature  dependence  of 
the  activation  energy  is  accounted  for  in  the  value  of  .  Thus,  there 
should  exist  a  value  Q,  for  a  specific  T0  which  is  constant  with  tempera¬ 
ture.  In  order  to  find  that  value  of  Q  an  iterative  procedure  was  used 
on  equation  (3.41)  for  a  trial  value  ofT0  .  Q  was  calculated  for  the 
different  temperatures  and  the  resulting  deviation  noted.  By  incrementing 
the  value  of  ^  ,  a  complete  set  of  deviations  was  obtained.  The  value 
of  Q  finally  chosen  corresponded  to  a  value  of  T0  which  gave  a  minimum 
variation  of  Q  with  temperature.  The  procedure  was  repeated  for  various 
moisture  contents  and  the  results  plotted  in  figure  3.8. 

The  energy  of  activation  required  to  rotate  a  water  molecule 

absorbed  on  0%  moisture  content  wood  is  found  to  be  about  13.4  to  14 

(51) 

Kcal/mole.  According  to  measurements  on  pure  ice  and  water  the 

activation  energy  for  ice  (3  bonded  water  molecules)  is  13.5  Kcal/mole 
and  for  water  the  bond  energy  is  4.5  Kcal/mole.  It  is  clear  from  figure 
3.8  that  the  activation  energy  required  for  polarization  of  pinewood 
drops  considerably  with  increasing  moisture  content.  However,  it  is  also 
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Fig.  3.8  -  Activation  energy  versus  moisture  content  for  Pine. 
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of  interest  to  note  that  close  to  fiber  saturation  (~30%M.C.)  the 
activation  energy  is  still  much  larger  than  that  required  to  break  one 
hydrogen  bond. 

The  factor  ^  is  usually  smaller  than  the  time  required  for  a  molecule 
to  go  from  one  equilibrium  position  to  another  because  of  internal 
excitation  of  the  molecule.  This  constant  further  includes  any  tempera¬ 
ture  dependence  of  the  activation  energy  as  pointed  out  previously. 

For  purposes  of  analysis,  T0  will  be  assumed  to  be  a  constant  with 
regard  to  temperature  and  moisture  content^ which  is  a  reasonable  approx- 
imation.  Hence  from  the  data  available  on  pinewood  a  value  of  T0  “  10" 
is  found.  Experimental  data  on  Douglas  Fir  and  Western  Hemlock  (presented 
in  chapter  5)  confirm  the  values  of  the  constants  found  here  from  pinewood 
data. 

Unfortunately  no  closed  form  solution  is  known  for  equations  (3.25), 
(3.27),  (3.28)  and  (3.29)  and  thus  numerical  integration  had  to  be  used 
to  solve  the  problem.  Using  I.B.M.'s  direct  access  APL  language  the 
problem  was  programmed  and  solved  to  a  high  degree  of  accuracy  on  an  I.B.M. 
360-67  digital  computer.  The  density  of  distribution  denoted  by  b  was 
obtained  by  iteratively  fitting  the  theoretical  curves  to  the  experi¬ 
mental  curves  to  obtain  a  fit  to  within  10%  accuracy.  Hence  theoretical 
plots  of  dielectric  constant  vs.  temperature  and  frequency  were  obtained 
for  semi-empirical  values  of  the  constants  b,  T0  »  Q>  and  ta  . 

Generally  the  agreement  between  theoretical  and  experimental  results 


was  very  satisfactory. 


. 


CHAPTER  4 


MEASUREMENT  METHODS  ANT)  TECHNIQUES 


INTRODUCTION 


When  looking  for  an  appropriate  technique  for  measuring  the  complex 
permittivity  of  a  material,  one  discovers  that  it  is  not  an  easy  task. 


At  microwave  frequencies  alone,  more  than  30  different  measuring  methods 
.  (23,  31,  32,  34) 


exist 


and,  depending  on  the  particular  material  to  be 


measured,  one  of  the  existing  methods  has  to  be  chosen,  after  which  the 
actual  exper imental  technique  is  often  modified  to  suit  the  particular 
need.  Some  important  factors  which  wiJ.l  affect  the  choice  are  a) 
accuracy  desired;  b)  simplicity  of  the  measurement;  c.)  simplicity  of  the 
calculations;  d)  geometry  and  make-up  of  the  material  to  be  tested; 
e)  environmental  requirements  for  the  sample;  f)  the  number  of  measurements 
to  be  made.  The  third  factor,  simplicity  of  calculations,  is  of  de¬ 
creasing  importance  with  availability  of  third  generation  computer 


facilities . 


In  the  present  investigation  it  was  necessary  to  vary  the  sample 
temperature  over  a  considerable  range,  which  was  one  of  the  reasons  for 
choosing  a  short-circuit  technique.  Short-circuited  mounts  have  only 
one  open  port  and  consequently  the  temperature  distribution  can  be  made 
more  uniform  than  far  a  two  port  device.  Moreover,  the  frequency  of 
interest,  namely  2450  MHz,  and  the  structure  of  the  material  to  be 
measured  (wood),  made  the  waveguide  technique  preferential.  The  fact 
that  nearly  1500  measurements  had  to  be  made  narrowed  the  choice  of  methods 
down  even  further  to  those  techniques  with  the  simplest  experimental  pro- 
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cedures.  Left  were  either  Von  Hippel’s  short-circuit  technique^"^  f 
figure  4,1,  or  the  bridge  method  developed  by  Buchanan  see  figure  4,2, 
Because  of  the  much  simpler  calculations  and  ease  of  measurement  in  the 
latter  method ,  it  was  the  measuring  technique  chosen,  /mother  important 
advantage  of  the  Buchanan  technique  is  that  front  face  reflections  of 
the  sample  do  not  affect  the  experimental  accuracy  as  they  do  in  Von 
Hippel’s  method, 

BASIC  METHOD 

Consider  the  sample  mount  shown  in  figure  4,3,  The  sample  is  backed 
by  a  short-circuit  plate  and  two  ports  in  the  broad  wall  of  the  waveguide 
are  designated  by  A  and  B,  The  distance  from  the  short-circuit  plate  to 
B  is  exactly  twice  the  distance  to  A, 

For  propagation  of  the  fundamental  mode  in  a  dielectric  filled  guide 


the  constant,  J7,  is  given  by_(31  32). 

Y  =  '  i71/  >■«}  V £  “P 

With  £,  the  complex  permittivity,  defined  by 

(4.1) 

£r£'+-j£"  and  P  =  (  /Xc}2 

(4.2) 

from  which 

£‘=p  +(Xo/2TTj2fP1-«V 

(4 .3) 

£"=  l(),0/zn)zc<s 

(4.4) 

where  A0  is  the  free  space  wavelength,  J3  and  are  the  phase  and  atten¬ 
uation  constants  of  the  material  and  A,  is  the  cut-off  wavelength  for  the 
waveguide.  In  order  to  find  values  for  £  and  ,  expressions  for  ^ 
and  J5  in  terms  of  the  experimental  quantities  need  to  be  obtained. 
Writing  the  expressions  for  the  electric  field  strength  in  the  material 
at  points  A  and  B  in  figure  4,3,  yields 


H 
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Fig.  4.3  -  Sketch  of  a  sample  mount 
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2^1  jcoi 
e 


The  ratio  of  equation  (4.6)  to  equation  (4.5)  gives 


(4.5) 


(4.6) 


JEi  =  (  Sinh  2X/ysinh  X/ -  z  cosh  H  (4.7) 

Experimentally  one  can  measure  Eg/E^  as  a  quantity  having  amplitude, r 
and  phase,©,  hence,  Equation  (4.8). 


H  _  2  cosh  W  = 

Solving  for  results  in 


V=  °c+jJ3  =  cash*'[(  r/z)ei91j 


(4.8) 


(4.9) 


Expanding  the  inverse  hyperbolic  cosine  function,  equation  (.4.9)  yields 


+  jJ3  =  Ln 


,  i©/2 

reJ  + 


-  (reje/Y 


(4.10) 


Hence,  equation  (4.10)  gives  the  value  of  oC  and  J3  in  terms  of  the  measured 
amplitude  and  phase  rand  9.  Because  equation  (4.10)  gives  a  multivalued 
solution,  it  is  necessary  to  repeat  the  measurement  for  a  different 
value  of  probe  spacing.  From  figure  4.3  it  is  clear  that  a  number  of 
ports  could  be  built  in  the  sample  mount  to  achieve  the  different  probe 
spacings.  However,  since  this  doubles  the  number  of  measurements  to  be 
taken  and  adds  experimental  complexity,  another  scheme  was  devised.  By 
attaching  the  sample  to  the  short-circuit  plunger,  the  distance  between 
the  short-circuit  and  the  pick-up  probe  could  be  set  to  any  desired 
value.  To  detect  the  position  of  the  probe  with  respect  to  the  short- 
circuit,  the  reference  arm  of  the  bridge  network,  see  figure  4.2,  was 
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shorted  out  by  means  of  a  short  circuit  switch.  The  short  circuit 
plunger  with  attached  sample  was- then  moved  all  the  way  in  past  the 
pick-up  probe.  By  watching  the  meter  indication  on  the  heterodyne 
detector  while  slowly  moving  the  plunger  back  out  of  the  guide,  the 
actual  short-circuit  position,  as  well  as  the  positions  of  the  minima 
and  maxima  of  the  standing  waves,  were  easily  determined. 

Since  the  short  circuit  position  is  moved  a  distance  /  away  from  the 
generator  for  the  second  measurement,  the  incident  wave  entering  the 
sample  will  have  advanced  its  phase  by  an  amount  equal  to 

(j)  — 

where  \c.  is  the  guide  wavelength  in  the  empty  guide.  Therefore,  the 
total  phase  difference  between  the  two  probe  positions  on  the  sample  is 
the  sum  of  the  phase  shifter  reading  plus  the  angle  0  given  above. 


If  the  principal  value  of  the  phase  angle  in  the  solution  of  Equation 
4.10  is  desired,  it  implies  that  the  length  -t  should  be  less  than  or 
equal  to  a  1/4  wavelength  in  the  material,  corresponding  to  the  first 
maxima  in  the  standing  wave  pattern  from  the  short-circuit  side.  The 
second  measurement  is  then  taken  at  twice  that  distance  or  at  the  first 
minima  from  the  short-circuit  side.  Ambiguity  in  the  measurement  is 
removed  by  the  above  procedure. 


Highest  accuracy  in  the  measurement  is  obtained  when  the  amplitude 
ratio,  r,  is  a  minimum  because  the  derivatives  of  Equation  (4.10)  taken 
with  respect  to  r  and  Q  are  directly  proportional  to  r.  By  means  of  the 
above  mentioned  procedure,  locating  the  minimum  of  the  standing  wave 
pattern,  corresponding  to  minimum  r,  was  a  simple  matter.  Therefore,  us* 
of  the  moving  short-circuit,  to  locate  the  1/4  and  1/2  wavelength 
positions  in  the  sample,  removes  the  multivalue  uncertainty  in  the 
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solution  of  equation  (4.10)  and  guarantees  maximum  accuracy. 

EXPERIMENTAL  APPARATUS 

t 

A  block  diagram  of  the  experimental  apparatus  is  given  in  figure  4.4. 

A  signal  from  the  microwave  generator  goes  through  a  directional  coupler 
where  part  of  the  signal  is  coupled  to  a  reference  arm.  The  reference 
arm  of  the  network  contains  a  precision  rotary  vane  attenuator  and  a  20 
db  variable  attenuator.  The  latter  is  used  for  setting  the  bridge  up 
for  maximum  sensitivity.  The  signal  in  the  main  arm  of  the  bridge  goes 

to  a  sample  mount  backed  by  a  moveable  short-circuit.  A  probe  in  the 

\ 

sample  mount,  flush  with  the  inside  of  the  guide,  picks  up  a  signal  of 
some  amplitude  and  phase.  The  unknown  signal  then  goes  through  a  precision 
phase  shifter,  is  mixed  together  with  the  reference  and  local  oscillator 
signals  and  the  output  fed  to  the  heterodyne  detector.  Both  attenuator 
and  phase-shifter  are  adjusted  to  obtain  a  deep  null  on  the  detector 
indicating  a  balanced  condition. 

To  reduce  the  number  of  experimental  runs,  the  signal  arm  of  the 
bridge  is  split  into  three  by  means  of  a  waveguide  junction,  each  arm 
of  which  goes  through  a  10  db  padding  attenuation , for  reducing  mount 
interaction, to  a  sample  mount  backed  by  a  moveable  short-circuit.  See 
Figure  4.4.  Only  one  of  the  sample  mounts  is  connected  in  the  bridge 
network  at  any  one  time  through  the  coaxial  selector  switch.  Precision 
dial-type  indicators  were  used  to  determine  the  short-circuit  position 
to  an  accuracy  of  0.001”.  Care  was  taken  to  ensure  that  both  arms  of 
the  bridge  were  properly  matched.  This  was  facilitated  by  using  a  swept 
frequency  ref . lactometer  and  thus  the  whole  system  was  aligned  for  a 
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Fig.  4.4  -  Experimental  apparatus 


62 


standing  wave  ratio  of  less  than  1.15.  Some  of  the  system  components 
will  be  discussed  in  more  detail  in  following  sections. 

PRECISION  PHASE  SHIFTER 


Since  no  satisfactory  phase  shifter,  coaxial  or  waveguide,  could  be 
obtained  commercially  for  2450  MHz,  it  was  necessary  to  design  a  phase 
shifter  particularly  suited  for  this  experiment.  The  specifications  were: 
a  minimum  accuracy  of  ±4°,  a  maximum  variation  in  insertion  loss  of  0.05 
db  over  a  range  of  0  to  360°,  and  a  direct  linear  readout  in  degrees. 
Specifications  were  met  by  utilizing  a  movable  short  circuit  plunger  in 
a  section  of  S-band  waveguide,  see  figure  4.5.  The  short  circuited 
guide  was  fed  through  the  back  of  the  short  by  means  of  a  coaxial  line 
and  coupling  loop  which  was  empirically  designed  to  achieve  a  good 
match  at  the  design  frequency.  However,  movement  of  the  short  circuit 
also  caused  the  coaxial  line  to  move  and  flex  which  resulted  in  excessive 

(33) 

spurious  phase  shift  even  though  Buchanan  ,  found  that  spurious  phase 
shift  with  his  equipment  was  negligible.  To  eliminate  this  problem,  a 
sliding,  constant  impedance,  coaxial  line  was  used  with  the  sliding 
section  attached  to  the  short  circuit  plunger.  If  now  the  waveguide 
short  circuit  position  is  changed,  the  length  of  both  the  coaxial  line 
and  of  the  waveguide  changes. 

Mathematically,  the  phase  shift,  ©3  ,  due  to  the  change  in  waveguide 
length  is  given  by 

©3  =  360  (  x/P\j)  (4’n)  _ 

where  X  is  the  distance  travelled  from  some  reference  point,  and  ^  is 
the  guide  wavelength.  The  phase  shift,  ©c  ,  due  to  the  change  in  coaxial 


.  ' 


<  - 


‘ 


6} 


64 


line  length,  is  similarly  given  by 


=  360  ( X/A0) 


(4.12) 


where  is  the  free  space  wavelength.  Since  the  travel  is  always  in 
the  same  direction  for  both  the  coaxial  line  and  for  the  waveguide,  one 
can  take  the  difference  between  equation  (4.11)  and  equation  (4.12)  to 
obtain  the  overall  phase  shift, 


(4.13) 


The  wavelength  in  the  waveguide  and  the  coaxial  line  is  fixed  for  any 
given  measurement  frequency  and  thus  the  total  phase  shift, 0  ,  becomes 
directly  proportional  to  the  position  of  the  short  circuit  plunger,  X  . 

The  total  length  of  travel  can  now  be  easily  calculated  by  setting  e  =  660 


equation  (4.13).  At  2450  MHz  in  S-band  guide,  WG284,  the  required  total 
travel  is  found  to  be  26.02  cm.  A  linearly  engraved  scale  was  made  to 
give  direct  readout  in  degrees  for  the  particular  test  frequency. 

It  is  important  to  insure  a  good  match  between  the  coaxial  line  and 
the  waveguide  in  order  to  keep  the  insertion  loss  variation  to  a  minimum. 
Note  also  that  the  phase  shifter  only  behaves  in  a  linear  fashion  if  the 
device  sees  a  matched  load  at  the  waveguide  end.  This  condition  is  quite 
easily  obtained  at  a  single  frequency.  After  fabrication  of  the  device, 
it  was  checked  for  linearity  and  accuracy  by  using  a  substitution  techniqu 
The  phase  shift  of  various,  accurately  known,  lengths  of  standard  guide 
were  measured  by  setting  the  phase  shifter  to  obtain  a  null  at  the  output 
with  and  without  the  waveguide  sections  in  the  circuit.  Results  indicated 
a  maximum  error  of  ±  1.5  degrees.  Maximum  variation  of  insertion  loss  of 


. 


'  •  ' 

. 


- 


65 


the  phase  shifter  was  measured  and  found  to  be  0.025  db.  Therefore,  this 
phase  shifter  design  was  subsequently  used  in  the  experiments. 

TEMPERATURE  CONTROL  OF  SAMPLE  MOUNTS 

To  obtain  information  on  the  temperature  behavior  of  the  dielectric 
properties  of  wood,  some  means  of  controlling  the  temperature  of  the 
sample  mounts  was  necessary.  The  low  temperature  limit  was  chosen  to 
be  -10°C  so  that  the  effect  of  freezing  temperatures  on  the  dielectric 
behavior  could  be  studied.  At  the  high  temperature  end  the  limiting 
factor  was  that  of  wood  decomposition,  hence  the  limit  was  specified  as 
100°C.  In  practice,  the  boiling  action  of  water  in  the  wood  cells 
dictated  a  maximum  attainable  temperature,  at  constant  moisture  content, 
of  about  95°C  at  3000  ft.  elevation.  Temperature  control  was  specified 
to  be  ±  1/2°C  to  give  acceptable  accuracy  on  the  temperature  behavior 
of  the  materials  under  test.  A  temperature  control  system  was  designed 
to  meet  the  specifications  and  is  described  below.  The  block  diagram  of 
the  apparatus  is  shown  in  Figure  4.6. 

Block  1  is  a  3/4  HP  refrigeration  unit  used  to  provide  cooling  of 
the  thermostatic  fluid.  It  was  capable  of  pulling  the  whole  system  down 
to  -10°C  when  all  components  were  well  insulated.  A  solenoid  in  the  re¬ 
frigerant  line  controled  the  action  of  this  equipment.  The  solenoid,  in 
turn,  was  governed  by  a  solid  state  control  unit,  which  received  its 
control  signals  from  a  Copper— Constant an  thermocouple  immersed  in  the 
thermostatic  fluid.  The  refrigerant  passed  through  a  50  ft.  copper  cool¬ 
ing  coil,  block  2,  completely  immersed  in  the  thermostatic  fluid,  thus 
ensuring  good  thermal  contact  between  the  two  media.  Heat  was  supplied 
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by  two  750  watt  immersion  type  heaters,  block  3,  which  also  were  con¬ 
trolled  by  the  control  unit  in  block  6.  The  sensing  tip  of  the  ther¬ 
mocouple  was  immersed  in  the  tank  containing  the  heaters  in  close 
proximity  to  the  intake  of  the  sample  mounts.  The  tank  was  wrapped  on 
the  outside  with  two  500  watt  sections  of  heating  tape  which  could  be 
turned  on  or  off  manually  as  required  to  maintain  a  given  heat  input  to 
the  system.  The  voltage  on  these  heating  tapes  could  be  controlled  by 
means  of  a  variac  to  obtain  some  control  over  the  total  heat  generated 
by  these  units. 

Ethylene  glycol  (antifreeze  without  additives)  mixed  with  water  was 
used  as  the  thermostatic  fluid.  It  has  good  thermal  properties,  is 
relatively  inexpensive,  and  the  mixture  has  a  boiling  point  of  about 
150°C.  Oil  was  tried  but,  because  of  the  high  circulation  rate  required 
to  maintain  accurate  temperature  control,  it  caused  excessive  foaming 
and  air  trapping.  The  system  was  operated  at  atmospheric  pressure  ana, 
therefore,  foaming  took  place  very  easily.  Block  4  represents  a  1  KP 
pumping  unit  capable  of  pumping  against  a  40  ft.  head  of  water  at  4.5 
gpm.  This  capacity  was  necessary  in  order  to  offset  the  losses  and 
accumulation  of  head  in  the  hoses,  valves,  etc. 

Three  identical  sample  mounts,  block  5,  were  fed  simultaneously  by 
means  of  a  set  of  manifolds,  each  of  which  had  twelve  ports.  The 
sample  mounts  consisted  of  8M  sections  of  waveguide  with  four  separately 
fed  water  jackets  soldered  to  each  side  of  the  guide.  The  water  jackets 
were  made  up  from  a  series  of  1/4"  copper  tubes  laid  parallel  to  each 
other  and  joined  by  manifolds  at  each  end.  In  this  way,  uniform  tempera¬ 
ture  distribution  across  the  mount  was  obtained.  Sample  mount  temperature 
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Fig.  4.7  -  Control  Circuit  for  Maintaining  Constant  Sample  Mount 
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was  continually  monitored  by  means  of  a  thermocouple  junction  soldered  to 
the  outside  waveguide  wall. 

The  solid  state  control  unit  in  block  6  consisted  of  a  series  of 
operational  amplifiers  and  Triacs,  three-terminal  solid  state  ac  switches, 
which  did  the  switching  and  controlling  of  the  heating  and  cooling  units. 

A  detailed  schematic  of  the  electronics  is  given  in  Figure  4.7. 

SAMPLE  MOUNT  DESCRIPTION 

A  sketch  of  the  actual  sample  mount  with  a  variable  short  circuit  is 
shown  in  Figure  4.8.  Mounted  in  the  center  of  the  broad  wall  and  flush 
with  the  inside  of  the  guide  is  a  lengthened  BNC  connector  which  is  used 
as  the  pick  up  probe.  Its  presence  in  the  waveguide  wall  caused  negligible 
disturbance  to  the  fields  in  the  guide.  At  the  open  end  of  the  guide,  a 
thin  mica  window  was  provided  to  minimize  moisture  loss  from  the  sample 
mount,  especially  at  the  higher  temperatures.  The  wTater  jacket’s  con¬ 
struction  is  evident  and  was  described  previously. 
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Fig.  4.8  -  Sample  Mount  Sketch  Showing  Cooling  Jacket  and  Probe. 
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To  obtain  a  good  variable  short  circuit,  a  resonant  box- type  short 
was  used  which  had  its  electrical  short  circuit  position  very  close  to 
the  physical  front  face  position  of  the  short.  Design  of  this  short  is 
described  in  reference  (35)  and  can  be  obtained  commercially  as  well  . 
To  accurately  position  the  short  circuit,  a  large  dial  type  position 
indicator  was  attached  to  the  back  of  the  short  making  it  possible  to 
read  the  short  circuit  position  to  0.001"  and  had  the  added  advantage 
of  easy  readability. 

CONDITIONING  OF  SAMPLES 


Moisture  content  in  wood  is  dependent  upon  temperature  and  humidity, 
hence,  to  condition  the  wood  specimens  it  is  necessary  to  control  these 
two  variables  in  some  way.  Furthermore  it  is  important  to  note  that  in 
order  to  get  a  uniform  moisture  content  in  a  specimen  a  fairly  long 
conditioning  period  is  necessary,  since  diffusion  of  the  moisture  between 
a  controlled  atmosphere  and  the  specimen’s  interior  is  a  slow  process. 

To  obtain  the  desired  data  for  a  plot  of  dielectric  constant  versus 
moisture  content  it  was  decided  to  use  specimens  of  the  same  wood  species 
conditioned  at  eight  different  equilibrium  moisture  contents,  covering 
the  range  of  0-25%. 

Moreover,  because  of  the  anisotropic  nature  of  wood,  due  to  its 
grain  structure,  it  would  be  useful  to  measure  the  dielectric  properties 
for  all  three  grain  directions,  namely,  tangential,  radial  and  longitudinal. 
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Because  wood  samples  are  non-uniform  in  structure  and  make-up,  three 
samples,  of  the  same  density  and  matched  as  close  as  possible  in  grain 
structure,  were  to  be  measured  under  identical  conditions  and  the 

I 

results  averaged.  Thus  it  follows  that  to  obtain  a  meaningful  set  of 
measurements  for  each  species  it  would  be  necessary  to  have  a  total  of 
8  (no.  of  different  moisture  contents)  times  3  (no.  of  grain  directions) 
times  3  (no.  of  samples  for  averaging),  or  72  samples  made  up.  It  was 
decided  to  double  this  number  to  allow  for  spoiled  samples,  initial  test 
runs,  and  some  added  tests  on  the  samples  if  desired. 

Two  species,  Western  Hemlock  and  Douglas  Fir,  were  chosen  for  the 
measurement  to  compare  the  effect  of  the  species  on  the  dielectric 
properties.  This  brought  the  total  number  of  samples  to  be  prepared  to 
288. 

In  order  to  be  able  to  condition  all  of  these  simultaneously  at  the 
required  moisture  contents  it  became  necessary  to  build  a  set  of  cabinets 
with  a  controlled  air  atmosphere.  Figure  4.9  shows  a  system  layout  for 
one  of  the  conditioning  cabinets.  Seven  of  these  units  were  placed  side 
by  side, each  with  its  own  control  circuit  for  humidification.  The  de- 
humidifying  part  of  the  system  feeds  the  seven  cabinets  in  series  with 
the  option  of  manually  controlling  the  dry  air  intake  in  each  cabinet  by 
means  of  a  valve.  All  cabinets  were  operated  at  their  own  ambient 
temperature  which  was  usually  somewhat  above  the  room  temperature  because 
of  heat  added  to  the  cabinet  by  the  steam  injection  process. 

When  the  humidity  in  the  cabinet  dropped  below  its  set  point,  the 
resistance  of  the  sensor  increased  and  caused  the  amplifier  to  put  out  a 
trigger  signal  causing  the  solid  state  relay  (Triac)  to  close  allowing 
the  vaporizer  to  generate  steam.  The  steam  entering  the  cabinet  was 
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Humidity  controlled  cabinets  -  System  layout 
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circulated  by  a  small  electric  fan  which  ran  continuously.  As  soon  as 
the  sensor  reached  the  required  equilibrium  value,  the  output  of  the 
differential  amplifier  became  zero  and  the  relay  opened,  thus  cutting 
off  the  power  to  the  vaporizer.  The  dehumidifier  ran  continuously,  but 
was  necessary  only  for  those  cabinets  which  were  set  for  humidities 
close  to  or  below  that  of  the  room  humidity.  The  door  seals  on  the 
higher  humidity  cabinets  were  purposely  made  leaky  so  that  overshoot 
effects  of  the  vaporizers  were  minimized.  An  on-offratio  of  1/5  to 
1/2  was  necessary  to  achieve  accurate  humidity  control,  depending  on  the 
effectiveness  of  the  dehumidifier  unit.  The  dehumidifier  sucked  moist 
air  from  the  cabinets  over  a  set  of  cooling  coils  which  were  kept  close 
to  0°C,  whereby  the  moisture  in  the  air  condensed.  The  dried  air  was 
passed  over  a  small  resistance  heater  raising  the  air  temperature  to 
that  of  the  cabinets. 

The  above  scheme  was  found  to  be  very  effective  in  controlling 
humidities  between  10%  and  90%  R.H.  to  an  accuracy  of  about  ±^%.  The 
Honeywell  humidity  sensors  do  tend  to  drift  somewhat  with  age,  but  a 
periodic  change  of  the  calibration  resistors  in  the  differential  ampli¬ 
fiers  would  bring  the  units  back  to  their  original  operating  points. 

This  procedure  was  found  to  be  necessary  about  once  every  six  weeks. 

This  time  interval  is  normally  sufficient  to  bring  a  reasonable  sized 
wood  sample  to  its  equilibrium  moisture  content. 

Figure  4.11  shows  a  schematic  of  the  solid  state  relay  used  for 
control  purposes.  These  units  have  proven  highly  satisfactory  in  a 
continuous  operating  environment.  Over  a  year  of  service  these  units 
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Notes 

14  Inverter  runs  at  about  3 -KHz ,  square  wave . 

2,  Trigger  requirement:  +12V  at  20mA;  minimum^ 8 ,57  Du, 
maximum- 24V  DC, 


3,  Tl-Hammond  Transf ormer ,  144G 
T2-Hammond  Transformer,  1441 
TRl-G.E.  Triac,  SC4 6B 

Fi rr  A  11  —  Solid  State  Relay  7  ,5 A,  117  V,AC, 
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have  cycled  many  millions  of  times  and  never  failed.  Mechanical  relays 
were  tried  initially,  but  were  found  highly  unreliable  and  caused  many  a 
system  failure  due  to  the  arcing  and  burning  of  contacts.  The  solid 
state  units,  on  the  other  hand,  have  proven  reliable,  relatively  inexpen¬ 
sive,  multi-purpose,  as  well  as  being  noiseless. 

SAMPLE  FITTING 

In  order  to  fit  the  samples  into  the  waveguide  the  raw  samples  had 
to  be  planed  or  milled  down  to  slightly  less  than  the  internal  waveguide 
dimensions.  The  raw  samples  which  were  about  30%  oversize  were  planed 
down  to  size  after  they  had  been  properly  conditioned.  Final  dimensions 
of  the  samples  were  nominally  2.825"  by  1.325"  by  3".  The  blocks  were 
then  returned  to  the  conditioning  cabinets  for  final  conditioning.  This 
left  a  clearance  of  .015"  between  sample  and  waveguide  walls  to  allow 
for  movement  of  the  samples  in  the  guide  ever  the  complete  temperature 
range.  The  actual  clearance  of  each  block  in  its  mount  was  measured 
just  prior  to  the  test  and  corresponding  corrections  were  applied  in  the 
final  calculations  of  the  dielectric  properties. 

* 

Double  sided  adhesive  transfer  tape  was  used  to  attach  the  samples 
to  the  short  circuit  plungers.  A  layer  of  tape  was  applied  to  the 
plunger  and  to  the  sample  end  face  which  had  negligible  effect  on  the 
measurement.  By  applying  reasonable  hand  pressure  the  blocks  were 
forced  against  the  plunger  face,  resulting  in  a  satisfactory  bond. 


*  Scotch  Brand  Adhesive  Transfer  Tape,  No.  465 

Minnesota  Mining  and  Manufacturing  of  Canada  Ltd. 
London,  Ontario 
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A  description  of  a  typical  measurement  using  the  bridge  technique, 

see  figure  4.4,  is  given.  A  sample  was  weighed  immediately  after  it 

was  removed  from  its  conditioning  cabinet,  its  weight  recorded  after 

which  adhesive  tape  was  applied  to  one  end  of  the  block  and  to  the  front 

face  of  the  short  circuit  plunger.  The  clearance  between  the  sample 

and  the  broad  waveguide  i^alls  was  Measured  and  recorded.  The  sample 

was  pressed  against  the  short  circuit  plunger  and  the  whole  assembly 

fitted  onto  the  sample  mount  and  the  sample  pushed  up  tight  against  the 

mica  window  to  minimize  air  circulation  around  the  block.  After  setting 

the  temperature  control  to  the  desired  value  the  whole  system  was  allowed 

0 

to  reach  an  equilibrium  temperature.  Allowing  30  minutes  per  10  C 
change  in  temperature  was  found  sufficient  for  stabilization.  To 
determine  the  time  interval  required  to  obtain  equilibrium,  the  bridge 
was  balanced  corresponding  to  a  null  at  the  detector.  When  no  further 
change  in  the  null  condition  of  the.  bridge  cccured  the  system  was  assumed 
at  equilibrium. 

By  means  of  a  thermocouple  soldered  to  the  waveguide  wall  the  mount 
temperature  was  continually  monitored.  Because  of  heat  losses,  the 
sample  temperature  was  usually  somewhat  different  from  the  mount  tempera¬ 
ture.  A  curve  of  sample  vs.  mount  temperature  was  obtained,  see  figure 
4.12,  by  inserting  a  thin  thermocouple  through  the  waveguide  wall  into 
the  center  of  the  sample,  after  equilibrium  was  reached.  This  procedure 
was  repeated  for  various  temperatures  in  the  range  and  for  a  few  different 
samples.  The  resultant  curve  was  a  straight  line.  Sample  temperature 


and  mount  temperature  were,  of  course,  equal  at  room  temperature. 
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Fig.  4.12  -  Sample  Temperature  versus  Mount  Temperature 
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equation  of  this  correction  curve  was  determined  and  since  a  computer 
was  used  for  calculations  the  temperature  correction  was  applied  to  the 
raw  input  data  as  part  of  the  computer  program. 

Once  equilibrium  was  reached,  the  actual  microwave  measurement  could 
begin.  With  the  waveguide  switch  closed  the  short  circuit  position  was 
adjusted  to  obtain  the  standing  wave  minimum  closest  to  the  physical 
short  circuit  position.  The  detector  would  show  a  deep  null  at  that 
point.  The  corresponding  reading  on  the  dial  indicator  was  recorded  as 
say,  X0  after  which  the  short  circuit  was  moved  away  from  the  pickup 
probe  till  the  second  minimum  in  the  standing  wave  pattern  was  detected. 
The  second  indicator  reading  was  recorded  as  .  Opening  the  waveguide 
switch,  the  attenuator  and  phase  shifter  were  adjusted  to  obtain  a  deep 
null  on  the  detector  and  the  attenuation  and  phase  shift  values  recorded. 
The  position  of  the  sample  and  short  circuit  was  then  changed  to  a  point 
midway  between  the  values  of  and  and  the  measurement  repeated. 

Half  the  difference  between  and  Xc  corresponded  to  the  distance  t?  as 
described  in  the  theory.  After  a  measurement,  the  sample  was  again 
moved  right  up  against  the  mica  window  and  the  temperature  control  set 
to  a  new  temperature.  To  run  a  complete  test  from  -10  to  1C'0°C  took 
about  8  hours  allowing  for  a  larger  initial  stabilization  time  since  the 
samples  originally  were  at  about  25°C.  By  using  three  sample  mounts 
three  tests  could  be  run  in  almost  the  same  amount  of  time,  since  the 
actual  measurement  time  was  short  compared  to  the  30  minute  stabilizing 


period . 
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MEASUREMENT  ERRORS 


A.  Short  circuit  positioning 

For  a  waveguide  short  circuit  plunger,  the  electrical 
short  circuit  position  is  often  somewhat  different  from  the 
physical  short  circuit  position.  In  the  case  of  the  resonant 
box- type  short,  used  in  the  experiments  described,  the 
electrical  short  circuit  position  was  found  to  lie  about  .9 
mm  ahead  of  the  physical  short  circuit.  However,  because 
of  the  technique  used  in  the  measurement  this  discrepancy 
was  of  no  real  concern  in  that  the  actual  standing  wave 
minima  were  located.  Moreover,  even  for  position  errors  of 
1  mm  (=  .040")  the  error  in  the  values  of  the  calculated 
dielectric  constant  was  less  than  1%.  Since  the  dial 
indicators  could  be  read  to  .001",  the  errors  due  to  the 
uncertainty  in  the  short  circuit  position  were  neglected. 

B.  Sample  clearance 

When  the  waveguide  is  not  completely  filled  with  a 
dielectric  material,  the  measured  value  of  dielectric  con¬ 
stant  will  be  less  than  the  true  value.  If  the  sample 
height  is  less  than  the  narrow  waveguide  dimension,  b,  the 
true  value  of  the  complex  dielectric  constant  is  given  by 


<4  *14) 


■where  c  is  the  clearance,  . 

The  effect  of  clearance  between  the  sample  and  the  narrow 


waveguide  walls  is  negligible  for  reasonably  small  clearances. 
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The  clearance  error  was  corrected  for  in  the  calculations  by 
use  of  equation  (4.14). 

C .  Temperature  variation 

The  sample  temperature  could  be  held  to  within  ±%°C . 

And  the  resultant  error  in  the  dielectric  constant  was  found 
to  be  less  than  1%. 

D .  Moisture  content  variation 

Variation  in  sample  moisture  content  between  the  low 
and  high  temperature  end  of  the  measurement  range  was  in¬ 
evitable.  No  humidity  control  inside  the  sample  mount  could 
be  provided  and  hence  at  the  higher  sample  temperatures 
some  of  the  moisture  was  driven  from  the  sample.  The  moisture 
partly  condensed  on  the  mica  window',  partly  filled  the  air 
space  in  the  mount  with  vapor,  and  partly  escaped  through 
the  back  of  the  short  circuit.  The  condensation  problem  was 
almost  eliminated  by  providing  a  heat  source  on  the  other 
side  of  the  mica  window,  thus  reversing  the  temperature 
gradient.  Only  at  temperatures  above  70°C  did  the  moisture 
loss  introduce  significant  errors  which  were  estimated  to 
lie  between  ±1  and  ±5%. 

E .  Mismatch  errors 

In  the  experiment  standing  wave  ratios  were  minimized 
throughout  and  wTorst  SWR  was  less  than  1.20.  However,  mis¬ 
match  errors  in  the  unknown  arm  of  the  bridge  tend  to  cancel 
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since  the  ratio  of  electric  fields  is  measured  in  the  sample 
and  the  reflection  coefficients  during  a  measurement  remain 
constant.  Therefore  this  source  of  error  was  found  to  be 
negligible. 

In  Appendix  A  an  analysis  is  given  for  the  calculation  of  errors 
in  the  dielectric  constant  due  to  experimental  variations  in  the  magnitude 
and  phase  of  the  detected  signals.  Variations  are  caused  by  fluctuations 
in  temperature  and  moisture  content  of  the  samples  and  by  uncertainties  in 
the  phase  shifter; attenuator  and  dial  indicators. 

Resultant  worst-case  errors  are  summarized  in  Table  4.1.  These  errors 
were  based  on  a  value  of  £  =2  and  £"=2.  Most  of  the  experimental  data 
gave  larger  values  of  £ ’  and  smaller  values  of  £"  which  would  decrease  the  % 
errors  considerably. 

Table  4,1  -  Summary  of  experimental  errors 


Amplitude 
and  phase 

Amplitude 
and  phase 

error  using: 
AA=  ±.05db, 
AQ=  ±2° 

length 

measurement 

error 

A-i-t.ozs  cm 

temperature 

variation 

error 

AT=±  .5  °C 

moisture 
content  var¬ 
iation  error 
up  to  ~/o  eC 

Total 

measurement 

error 

A£' 

<*•!% 

<ti% 

<±2% 

<i  sK 

A  £" 

<t  i% 

<1  1.5% 

<  1 1  % 

<ts% 

A  tan  6 

<  -+  ^  % 

<ti% 

<±l  % 

<tl% 

During  the  course  of  this  work  a  new  Swept  Frequency  Dielectric  Measure- 

07) 

ment  Technique  was  developed  by  Tinga  and  Edwards  .  This  technique 
proved  to  be  fast  and  reasonably  accurate.  However,  samples  measured  by 
this  technique  had  to  be  a  few  wavelengths  long.  Since  the  samples  had 
already  been  cut  to  a  shorter  length  and  conditioned 5the  Swept  Frequency 
Method  could  not  be  used  efficiently  for  all  the  measurements,  but  it  was 
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used  to  check  some  of  the  results  on  Douglas  Fir  previously  measured  by 
the  bridge  method.  This  procedure  confirmed  the  results  obtained  by  the 
bridge  technique  to  within  ±5%. 

All  experimental  results  obtained  by  the  procedures  described  here  are 
presented  in  graphical  form  in  chapter  5. 


f 


CHAPTER  5 
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CALCULATIONS  AND  RESULTS  FOR  TEMPERATURE  DEPENDENCE 

Following  the  procedure  outlined  in  chapter  4  the  dielectric 
properties  of  Douglas  Fir  were  measured,  and  the  data  processed  on  a 
digital  computer*  Various  plots  of  the  important  variables  are  pre¬ 
sented  as  functions  of  temperature,  moisture  content  and  frequency. 
Theoretical  results  are  compared  with  experimental  data,  and  the  phy¬ 
sical  meaning  of  the  relaxation  process  is  discussed  more  fully. 

A.  Experimental  Results 

1*  Dielectric  Constant  Versus  Temperature 

In  figures  5.1  through  5.9  the  variation  of  the  complex 
dielectric  constant  with  temperature  is  given.  All  measurements 
were  taken  at  a  frequency  of  2450  MHz,  which  was  chosen  for  the 
simple  reason  that  it  is  the  centre  frequency  of  a  band  designated 
for  industrial,  medical  and  scientific  uses.  It  so  happens  that  at 
this  particular  frequency,  over  the  temperature  range  shown  the 
dielectric  loss  curves  for  wood,  exhibit  distinct  maxima.  Referring 
back  to  the  theory  of  chapter  3  this  type  of  behavior  is  predicted 
to  occur  quite  generally.  However,  since  some  of  these  data  were  ob¬ 
tained  before  chapter  3  was  developed  it  was  at  first  not  obvious  that 
these  maxima  would  occur  at  this  frequency  in  the  rather  narrow  temp¬ 
erature  range  of  0-100 °C.  For  this  reason  many  of  the  tests  were 
repeated  and  the  results  shown  to  be  consistent  with  the  previous 
ones.  The  results  are  very  revealing  about  the 
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Fig  5.6  -  Loss  tangent  of  Douglas  Fir  asa  function  of 

temperature  at  various  moisture  contents  measured  at  2450 
MHz.  with  E-field  in  the  tangential  grain  direction. 
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Fig.  5.8-  Dielectric  loss  of  Douglas  Fir  as  a  function  of 

temperature  at  various  moisture  contents  measured  at  2450 
MHz.  with  E-field  in  the  radial  grain  direction. 
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Fig.  5.9  -  Loss  Tangent  of  Douglas  Fir  as  a  Function  of  Temperature  at  Various  Moisture  Contents  Measured 
at  2450  MHz  with  E-field  in  the  Radial  Grain  Direction. 
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nature  of  wood  as  a  dielectric.  Maximum  loss  temperatures, 

Tm,  are  plotted  for  Douglas  Fir  in  figure  5.10  for  the  longitudinal 
grain  direction.  From  the  graph  it  is  apparent  that  the  maximum 
loss  temperatures  of  £"  and  tan  6 are  quite  different.  Theoreti¬ 
cally,  for  low  dispersion  dielectrics  (£0  -Co  small) ,  the  maxima 
f or  £ "  and  tan  6  should  occur  at  nearly  the  same  temperature  and 
frequency.  As  dispersion  and  loss  increase,  the  two  curves  diverge, 

as  is  evident  from  figure  5.10.  A  similar  plot  of  T  versus  moisture 

m 

* 

content  is  given  for  Western  Hemlock  in  figure  5.11  . 


Only  small  variations  in  density  were  found  in  the  samples 
used  for  the  tests  and  corrections  for  these  variations  were  small 
and  hence  neglected.  However,  the  correction  for  the  real  part  of 


the  dielectric  constant  is  linearly  proportional  to  density  and 
further  data  on  such  corrections  may  be  found  in  the  literature 


(50) 


_2 .  Calculations  of  Activation  Energies 

In  order  to  produce  theoretical  data  for  these  materials  it  is 
necessary  to  calculate  the  activation  energy  required  in  the  polar¬ 
ization  process.  Invoking  the  theory  of  chapter  3,  the  activation 
energies  are  calculated  by  the  use  of  equation  (3.32)  which,  when 
rewritten  becomes 

Q=  RTLn(l/wT0)  (5.1) 

T  is  assumed  to  have  the  same  value  for  all  wood  species  and  is 
'o 

—18 

taken  to  be  1  x  10  sec  (derived  from  Pinewood  data  in  chapter 

3).  To  evaluate  Q,  the  values  for  T  ,  T.  ,  k)  and  R  are  put  in 

m 

equation  (5.1).  Figure  5.12  shows  the  calculated  activation 
(53) 

*  A.  J.  Curtis  gives  a  plot  of  T  versus  percentage  composition 
for  polyvinylchloride-tricresylphosphate.  The  functional  relation¬ 
ship  there  is  found  to  be  similar  to  that  of  wood. 
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Fig,  5,10  -  Maximum  Dielectric  loss  temperature  as  a  function  of  moisture 
content  for  Douglas  Fir  measured  at  2450  MHa  -with  E-field  in 
the  longitudinal  grain  direction,  . 
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Fig,  5,11  -  Maximum  dielectric  loss  temperature  as  a  function  of  moisture 
content  for  Western  Hemlock  measured  at  2450  MHz  'with  E-field 
in  the  longitudinal  grain  direction', 
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energy  of  polarization  as  a  function  of  moisture  content  for  Douglas 
Fir,  Western  Hemlock  and,  for  comparison,  the  curve  for  Pinewood  is 
given  as  well. 


It  is  well  known  that  wood  is  structurally  an  anisotropic  sub¬ 
stance  which  results  in  an  anisotropic  behavior  of  its  mechanical 
and  electrical  properties.  It  turns  out  that  the  dielectric  constant 
for  wood  measured  with  the  electric  field  along  the  grain  (longi¬ 
tudinal)  is  considerably  higher  than  that  measured  in  the  radial 
and  tangential  directions.  From  the  data  obtained  this  can  be 
interpreted  as  partially  due  to  a  variation  of  activation  energy 
of  polarization  with  grain  direction.  Figure  5.12  shows  the 
activation  energy  for  Douglas  Fir  in  the  tangential  direction  to 
be  somewhat  higher  than  that  for  the  longitudinal  grain  direction. 

The  difference  in  slope  of  the  curves  in  figure  5.12  for  the 

different  wood  species  is  likely  related  to  the  crystallinity  of 
* 

the  material  ,  Western  Hemlock  and  Pine  both  are  less  crystalline 
than  Douglas  Fir  which  is  also  reflected  in  their  density  and 
hardness . 

3. .  Bound  Water  Concepts 

The  initial  activation  energy  involved  in  the  polarization 
process  of  wood  at  0%  moisture  content  is  found  to  lie  between  13.5 


An  increase  in  crystallinity  decreases  the. maximum  loss  and  shifts 
the  loss  maximum  to  higher  temperatures V  .  Decreased  crystallinity, 
hence  decreased  hydrogen  bonding,  results  in  an  increased  dielectric 
constant.  It  is  found  that  a  linear  relationship  exists  between 
dielectric  constant  ,^gd  g^e  ^gystallinity  index  for  various 
cellulose  materials v  *  * 


X 
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%  Moisture  Content 

Fig.  5.12  -  Activation  energy  of  polarization  versus  moisture  content  for 
Douglas  Fir,  Western  Hemlock  and  Pine  "with  E-field  in  the 
longitudinal,  tangential  and  radial  grain  directions. 
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to  14  Kcal/mole.  This  reveals  the  type  of  bonding  taking  place 

(51) 

between  the  cellulose  and  water  molecules.  Hasted  and  others 
(54) 

have  done  extensive  tests  on  the  dielectric  properties  of  ice 
and  water.  It  is  found  that  for  ice  the  activation  energy  of 
polarization  is  13.5  Kcal/mole  and  for  water  the  value  drops  to 
4.5  Kcal/mole.  Assuming  the  bond  bending  and  breaking  theory, 

(59) 

advanced  in  the  literature  ,  holds,  three  or  four  hydrogen 
bonds  occur  in  ice  and  for  water  only  one  single  hydrogen  bond  is 
present  between  molecules.  Hence,  to  break  a  hydrogen  bond  com¬ 
pletely  requires  about  4.5  Kcal/mole,  whereas  breaking  and  immediately 
reforming  of  a  bond  is  estimated  to  require  about  3  Kcal/mole. 
Moreover,  bond  bending  may  occur  in  which  case  the  required  energy 
lies  somewhere  between  0  and  3  Kcal/mole.  In  light  of  these 
facts,  it  is  reasonable  to  assume  that  the  initial  water  molecules 
absorbed  on  wood  are  held  in  a  fashion  very  similar  to  that  of  pure 
ice.  As  the  moisture  content  in  the  wood  increases  the  water 
molecules  are  bound  less  tightly  than  before  because  it  becomes 
possible  for  the  molecules  to  make  water  to  water  bonds  as  well, 
which  requires  a  lower  energy.  This  should  happen,  for  example, 
after  the  initial  monomolecular  layer  of  water  is  filled  up, 
which  occurs  at  approximately  6%  moisture  content 


Bound  water  phenomena  in  wood  have  been  discussed  by  Voss 


(57) 


from  measurements  at  3  cm  wavelengths  and  by  others^’  ^  . 


However,  actual  energies  associated  in  the  dielectric  polarization 
process  in  wood  were  never  given  explicitly.  From  the  data  and 
the  analysis  presented  here  it  is  now  possible  for  the  first  time 
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to  see  explicitly  the  relationship  between  intrinsic  material  parameters 

and  the  macroscopic  dielectric  properties  of  a  bound-water  system.  An 

extremely  important  observation  can  be  made  about  the  so  called  free _ 

water  in  a  substance.  It  has  often  been  assumed  that  water  absorbed 

beyond  a  certain  moisture  content,  for  example,  the  fiber  saturation 

point  in  wood,  behaves  not  as  bound  water,  but  as  ordinary  liquid 

water.  However,  Hasted presents  ample  data  to  show  the  relaxation 

process  of  liquid  water  requires  an  activation  energy  of  only  4.5 

Kcal/mole,  whereas  the  data  on  wood  indicate  that  at  fiber  saturation 

(26-30%  M.C.),  the  activation  energy  of  polarization  is  still  as  high 

as  10  Kcal/mole  (see  figure  >''5.12.)  Hence,  unless  a  sharp  break 

occurs  in  the  relationship  between  activation  energy  and  moisture  con- 
* 

tent  ,  it  is  concluded  that  bound  water  effects  extend  beyond  the 
fiber  saturation  point  of  wood.  Other  experimental  evidence  which 

(26) 

supports  such  a  behavior  is  found  in  a  paper  on  wool-water  systems v 


Detailed  arguments  as  to  the  bonding  mechanism  between  water, 

(59) 

organic  and  inorganic  molecules  can  be  found  in  the  literature 
The  discussion  given  there  allows  for  the  rather  high  activation 
energies  found  for  wood-water  systems.  Based  on  the  experimental  data 
given  in  this  thesis  it  can  be  assumed  that  the  relaxation  of  water 
in  high  moisture  content  wood  still  involves  2  or  3  bonded  water 
molecules,  and  hence  does  not  have  dielectric  properties  equal  to 
that  of  pure  liquid  water. 


All  experimental  evidence  points  to  a  smoothuf^ansition  of  functional 
relationships  at  the  fiber  saturation  point ^  .  Similar^arguments 

were  presented  in  a  quite  different  context  by  C.  Skaar 
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!L'  Dielectric  Constant  Versus  Moisture  Content 

Data  from  the  dielectric  constant  versus  temperature  curves  can 
be  used  to  obtain  plots  for  the  variation  of  dielectric  constant  with 
moisture  content  using  temperature  as  the  parameter.  This  is  done 
for  Douglas  Fir,  longitudinal  grain,  in  figures  5.13  >  5.14  ,  and 

5.15..  The  real  part  of  the  dielectric  constant  increases  monotoni- 
cally  with  moisture  content  up  to  at  least  the  fiber  saturation  point 
and  increases  with  temperature  as  well.  Dielectric  loss  on  the  other 
hand  shows  a  more  complex  behavior  which  is  due  mainly  to  the  tempera¬ 
ture  dependence  of  the  relaxation  times.  Thus  at  the  higher 
temperatures  the  loss  factor  tends  to  level  off  with  increasing  moisture 
content.  Only  a  few  curves  are  given  in  figure  5.14  for  the  sake  of 
clarity.  The  data  compare  favorably  with  some  data  obtained  by  James 

(34) 

and  Hamill  .  Their  results  were  obtained  at  "room  temperature" 
interpreted  here  as  23°C.  . 

From  the  experimental  data  on  Douglas  Fir  and  Western  Hemlock  it 
is  possible  to  find  the  relationship  between  the  values  of  the 
dielectric  constant  at  the  maximum  loss  temperatures  and  moisture 
content.  The  data  are  presented  in  Figure  5.16,  which  indicate  a 
linear  dependence  between  the  parameters.  These  curves  will  be  used 
to  calculate  the  constants  needed  for  some  theoretical  curves  presented 
in  the  next  section. 

B.  Theoretical  Results 

With  the  aid  of  the  experimental  data  presented  in  the  previous  section, 
the  equations  developed  in  chapter  3  can  be  solved  yielding  the  theoretical 
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Fig*  3*13  -  Dielectric  constant  of  Douglas  Fir  at  2430  MHz  as  a 
function  of  moisture  content  v/ith  temperature  as  the 
parameter  and  E-field  in  the  longitudinal  grain 
direction* 
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Moisture  Content  v/ith  Temperature  as  the  Parameter  with 
E-field  in  the  Longitudinal  Grain  Direction* 
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Fig ,  5 .15  -  Dielectric  Loss  Tangent  of  Douglas  Fir  at  2450  MHz  as  a 

Function  of  Moisture  Content  -with  Temperature  as  the  Parameter 
with  E-field  in  the  Longitudinal-  Grain  Direction. 
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i’ig,  5,16  -  Dielectric  Constant  and  Dielectric  Loss  at  the  Maximum  Loss 

Temperatures  as  a  Function  of  Moisture  Content  for  Douglas  Fir 
and  Western  Hemlock  Measured  at  2450  MHz  vdth  E~field  in  the 
Longitudinal  and  Tangential  Grain  Directions, 
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behavior  of  the  dielectric  constant  for  wood. 

Equations  (3.25)  and  (3.27)  were  solved  numerically  on  a  digital 

computer  for  various  moisture  contents.  This  required  values  of  Q,  and 

~18 

b.  The  values  used  were  T0  =  10  sec  ,  Q  as  found  from  figure  5.12 
and  b  was  given  some  trial  value.  The  low  and  high  frequency  dielectric 
constants  were  obtained  by  integration  of  equations  (3.33)  and  (3.36). 

The  resulting  solution  was  compared  •with  the  experimental  curves  for 
C(T^  and  b  was  incremented  until  a  fit  of  better  than  ±10%  between  the 
theoretical  and  experimental  dielectric  loss  was  achieved.  This  procedure 
was  carried  out  automatically  for  a  range  of  0-100  C.  Repeating  the 
calculations  for  each  moisture  content,  values  of  b  were  found  for  all 
the  curves.  It  turns  out  that  the  value  of  b  is  approximately  constant 
(b=.255)  with  moisture  content.  For  comparison,  the  experimental  and 
theoretical  values  of  £  ’  and  £ "  versus  T  for  a  moisture  content  of  17.5% 
are  plotted  in  figures  5.17  and  5,18.  At  the  higher  temperature  the 
experimental  values  for  the  real  part  of  the  dielectric  constant  are  too 
low  because  of  moisture  loss  during  the  measurement.  Since  £?  increases 
rapidly  with  moisture  content  and  temperature,  see  figure  5.13,  the  effect 
is  quite  noticeable.  However,  percentage  error  is  less  than  ±10%.  Curves 
for  the  other  moisture  content  values  give  similar  results,  but  are 
omitted  for  the  sake  of  clarity.  Because  the  theoretical  data  lie  so 
close  to  the  experimental,  any  further  calculations  will  be  done  on  the 
theoretical  data,  since  its  functional  relationships  are  known. 


Families  of  curves  for  the  theoretical  complex  dielectric  constant 
of  Douglas  Fir  as  a  function  of  temperature  are  presented  in  figures  5.19 
and  5.20,  with  moisture  content  given  as  a  parameter.  These  two  figures 
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Fig,  5,17  -  Theoretical,  and  experimental  plot  of  dielectric  loss  versus 
temperature  for  17.5$  m,c.  Douglas  Fir  -with  E-field  in  the 
longitudinal  grain  direction;  measured  at  2450  MHz, 


C 


r\ 


0 


— r~ 


i 

_  i  . 
10 


~io" .  50  ” . 40  ~~  ‘  50 

Temperature  °C 


k 


Theoretical 

Experimental 


60  ~  . 70  80  . SO 


Fig,  5,18  -  Theoretical  and  experimental  plot  of  dielectric  constant  versus  tern 
temperature  for  17,5$  m,c,  Douglas  Fir  vrith  E-field  in  the 
longitudinal  grain  direction;  measured  at  2450  MHz. 
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Fig.  5.19  -  Theoretical  dielectric  constant  versus  temperature  at  various 
moisture  contents  for  Douglas  Fir  at  2450  MHz  -with  ga-field 
in  the  longitudinal  grain  direction. 
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Fig.  5.20  -  Theoretical  dielectric  loss  versus  temperature  at  various 
moisture  contents  for  Douglas  Fir  at  2450  MHz  with  E-field 
in  the  longitudinal  grain  direction. 
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can  be  used  as  convenient  reference  charts  for  the  dielectric  constant  of 
Douglas  Fir.  Calculations  of  the  slope  of  the  above  two  graphs  yield  the 
temperature  coefficients  of  the  real  and  imaginary  parts  of  the  dielectric 
constant  .  Hence,  curves  of  d£?/dT  and  d£"/dT  versus  moisture  content  are 
plotted  in  figures  5.21  and  5.22  with  temperature  as  a  parameter. 


For  completeness,  the  theoretical  values  of  the  dielectric  constant 
versus  moisture  content  are  also  calculated  and  presented  in  the  form  of 
figures  5.23,  5.24  and  5.25,  with  temperature  as  the  parameter. 

All  of  the  above  results  are  valid  at  the  frequency  of  2450  MHz.  From 
the  same  experimental  data  however,  the  behavior  of  the  dielectric  constant 
at  other  frequencies  can  be  determined.  With  the  help  of  equation  (3.26) 
it  is  possible  to  predict  the  relaxation  times  of  wood  at  various  moisture 
contents,  thus 

Q/RT 

(3.26) 


r=r0e 

Table  5.1  -  Activation  energy  and  relaxation  time  for 
Douglas  Fir  versus  moisture  content. 


Moisture 
Content ,% 

Activation 

Energy,  Q 
Cal/moLe 

Relaxation  Time  in  sec,  (x  10*^)  at 
0°C  20°C  40°C  60°C  80°C 

: 

100°C 

.1 

13660 

8640 

1550 

346 

92.4 

28.7 

10.1 

.5 

13600 

7730 

1400 

314 

84.4 

26.4 

9.32 

1 

13500 

6420 

1180 

267 

72.6 

22.9 

8.14 

2 

13310 

4530 

849 

197 

54.5 

17.4 

6.3 

5 

12800 

1770 

354 

86.7 

25.2 

8.42 

3.17 

8 

12300 

704 

150 

38.8 

11.8 

4.13 

1.61 

10 

12000 

405 

89.5 

24.0 

7.52 

2.69 

1.08 

15 

11350 

122 

29.3 

8.43 

2.82 

1.07 

.448 

20 

10900 

53.3 

13.5 

4.09 

1.43 

.561 

.244 

25 

10550 

28.0 

7.41 

2.33 

.841 

.341 

.152 

30 

10270 

16.7 

4.58 

1.48 

.551 

.229 

.104 

In  the  process  of 

microwave 

heating,  the  temperature 

coefficients  are 

of  the  treated  dielectric. 
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Fig.  5.21  -  Temperature  coefficient  of  the  theoretical  dielectric 
constant  versus  moisture  content  at  various  temperatures 
for  Douglas  Fir  at  2450  MHz  with  E-field  in  longitudinal 
grain  direction. 
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Fig*  5*22  -  Temperature  coefficient  of  the  theoretical 
dielectric  loss  versus  moisture  content  at  various 
temperatures  for  Douglas  Fir  at  2^50  MHz  with  E-field 
in  longitudinal  grain  direction* 


Dielectric  Constant  £  x 


114 


Content  at  Different  Temperatures  at  2450  LIHz  for 
Douglas  Fir  mth  E-field  in  the  Longitudinal 
Grain  Direction. 
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Generally,  care  must  be  taken  in  the  calculation  of  relaxation  times  from 

(  fifTi 

temperature  varying  data  at  a  given  frequency.  Frohlicn  states  that 
(  C0  -  )  and  ( £  -  £o o)  must  be  known  as  functions  of  temperature  before  T 

can  be  calculated.  Now,  dE^/dT  is  proportional  to  the  linear  expansion 
coefficient  of  the  material,  see,  for  example,  chapter  3,  which  for  various 
wood  species  is  of  the  order  of  one  part  per  million.  As  a  result,  the 
temperature  dependence  of  the  high  frequency  dielectric  constant  is 
negligible.  Furthermore,  data  on  Douglas  Fir,  for  example,  figure  5.1, 
indicate  that  the  temperature  dependence  of  £0 ,  the  low  frequency  or, 
conversely,  the  high  temperature  value  of  £ ' ,  is  also  very  small  -  £' 
becomes  constant  as  T  increases.  It  is  concluded  that  values  of  T  may 
safely  be  calculated  by  the  direct  use  of  equation  (3.26).  When  previously 
defining  this  equation,  it  was  also  stated  that  both  T0  and  Q  were  assumed 
to  be  temperature  independent  which  according  to  the  above  paragraph  has 
experimental  justification. 

*“18 

Applying  equation  (3.26)  for  T0  =10  sec  and  for  values  of  Q  found 
from  figure  5.12,  the  relaxation  times  for  Douglas  Fir  are  evaluated  and 
summarized  in  Table  5.1.  Figure  5.26  is  a  graphical  presentation  of  the 
relaxation  time  of  Douglas  Fir  versus  moisture  content  with  temperature  as 
the  parameter.  This  log-log  plot  shows  quite  clearly  that  at  about  6% 
moisture  content  the  behavior  of  relaxation  time  with  increasing  moisture 
content  changes  to  a  nearly  linear  log-log  plot,  and  at  very  high  moisture 

t 

contents  (30-200%)  the  relaxation  time  approaches  that  of  pure  liquid 
water.  The  last  part  of  the  curves  is  shown  by  broken  lines  to  indicate  the 

I 

probable  asymptotic  approach  to  the  pure  water  relaxation  times  at  the  various 
temperatures^51^.  This  particular  graph,  figure  5. 26, is  very  useful  in 


•m 


. 


Relaxation  Time 
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probable  asymptotic  approach  to  the  pure  water  relaxation  time: 
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Vj 


Frequency  Hz 

Fig,  5#27  —  Theoretical  dielectric  constant  and  loss  versus  frequency  at  different  temperatures  and 

a  moisture  content  of  5%  for  Douglas  Fir  with  E-field  in  the  longitudinal  grain  direction 
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determining  the  temperature  behavior  of  wood  at  any  given  frequency  and 
moisture  content.  For  example,  if  it  is  desired  to  know  the  character¬ 
istics  of  Douglas  Fir  at  915  MHz  the  relaxation  time  corresponding  to  this 
frequency  is  l/2rrf  =  1.74  x  10  ^  sec.  Drawing  a  line  for  T  =  1.74  x  10 
parallel  to  the  moisture  content  axis  it  is  immediately  apparent  at  which 
temperatures  and  moisture  contents  the  maxima  in  the  dielectric  loss  curves 
occur.  In  the  example  cited  a  20%  M.  C.  Douglas  Fir  will  exhibit  a  loss 
maximum  at  about  18°C.  Use  of  this  figure  clearly  shows  that  to  shift 
the  loss  maxima  to  a  higher  temperature  for  a  fixed  moisture  content  the 
applied  frequency  must  increase.  To  obtain  positive  temperature  coefficients 
on  the  dielectric  constant  of  Douglas  Fir  from  0  to  30%  moisture  content  a 
frequency  of  8000  MHz  or  higher  needs  to  be  applied. 

Having  found  values  for  the  relaxation  time,  the  frequency  response  of 

the  dielectric  can  now  be  obtained  by  using  equations  (3.25)  and  (3.27)  for 

*  _ _ - 

various  values  of  frequency,  moisture  content  and  temperature  .  The  result¬ 
ing  frequency  varying  data  are  presented  in  figures  5.27  and  5.28  based  on 
the  activation  energy  of  polarization  for  longitudinal  Douglas  Fir. 

Figure  5.27  shows  the  dielectric  constant  as  a  function  of  frequency 
at  different  temperatures  for  longitudinal  Douglas  Fir  at  a  moisture 
content  of  5%.  In  figure  5.28  the  dependence  of  the  dielectric  constant 
on  frequency  at  20°C  for  various  moisture  contents  is  shown.  Qualitatively 
these  results  agree  very  well  with  data  on  Pinewood  published  by  Trapp 
et  al^^.  Their  dataiwere  measured  with  the  electric  field  in  the  radial 
grain  direction,  whereas  the  data  on  Douglas  Firwere  measured  in  bhe 


These  results  as  all  the  others  were  obtained  by  the  use  of  a  digital 
computer . 
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longitudinal  grain  direction.  Hence  the  magnitudes  of  the  dielectric 

constant  for  the  longitudinal  Douglas  Fir  are  higher  than  the  radial  Pine- 

wood.  At  the  low  frequency  end, conductivity  effects  and  charged  double 
(27) 

layers  tend  to  increase  the  dielectric  constant  considerably.  Since 

these  effects  were  not  accounted  for  in  the  treatment  given  here  the  values 
for  dielectric  constant  below  100  MHz  are  too  low.  Nevertheless  the  agree¬ 
ment  between  the  theoretical  curves  in  figures  5.26  and  5.27  and  experimental 
frequency  data  as  found  in  the  literature^)  confirms  the  basic  assumptions 
underlying  the  theory  given  in  this  thesis.  Furthermore,  one  can  make 
reasonable  predictions  of  the  dielectric  properties  at  other  frequencies 
in  the  range  between  100  MHz  and  20  GHz. 

By  the  use  of  equations  (3.35)  and  (3.39)  the  theoretical  values  of 

the  low  and  high  frequency  dielectric  constant  can  be  calculated  using  the 

empirical  values  f or  £ (T  )  as  given  in  figure  5.16.  The  resultant  dependence 

m 

of  E0  and  Eoo  on  moisture  content  is  given  by  figure  5.29. 

In  assessing  which  parameters  are  basic  to  obtain  the  dielectric 
constant  for  wood  it  follows  from  this  chapter  that  only  the  activation 
energy'  Q  as  a  function  of  moisture  content,  and  the  variation  of  the  low 
and  high  frequency  dielectric  constants  must  be  known  accurately.  From 
the  experimental  data  on  Douglas  Fir  the  following  empirical  relations 
for  activation  energy  of  polarization  was  determined  by  a  leash  squares  fit 
process .  (5.2) 

Q  =  2.716  M2-  194-7 M  +  13690 

where  M  is  thejmoisture  content  of  Douglas  Fir,  Equation  (5.2)  fitted  the 
data  to  better  than  1/2%. 


Low  and  High  Frequency  Dielectric  Constants 
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Similarly  for  the  low  and  high  frequency  dielectric  constant  the 
empirical  relations  were  found  to  be 


£0=  .1 9JJM  +  Z.4 


(5.3) 


-.0I63M+I.95* 


(5.4) 


Equations  (5.2),  (5,3),  and  (5.4)  are  valid  between  0  and  30%  moisture 
content . 

Conclusion 

Combination  of  equations  (3.25),  (3.27),  (5.2),  (5.3)  and  (5.4)  gives 
a  complete  description  of  the  dielectric  properties  of  Douglas  Fir  over  a 
wide  frequency  and  temperature  range  and  at  any  moisture  content  between  0 
and  30%.  This  defined  system  of  equations  though  tedious  and  difficult  to 
solve  by  hand  was  easily  solved  on  a  digital  computer  using  APL,  a  direct 
access  interpreter  type  language,  requiring  only  seconds  of  computer  time 
and  10  to  20  kilobytes  of  core  storage. 

The  mathematical  description  of  a  wood  system  has  the  great  advantage 


of  improving  predictability  of  a  microwave  wood  heating  system’s  performance 
under  operating  conditions  where  frequency,  temperature  and  moisture  content 
may  all  change  continually. 


■ 
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CHAPTER  6 

■APPLICATION  OF  MIXTURE  THEORY 

In  this  chapter  the  dielectric  mixture  relations  developed  in  chapter 
2  will  be  applied  to  a  cellulose-air-water  mixture  such  as  Douglas  Fir, 
Values  of  the  complex  dielectric  constant  will  be  calculated  as  a  function 
of  moisture  content  and  frequency  and  the  results  compared  with  experimental 
data. 

WOOD  AS  A  MIXTURE 

Wood  is  a  complex  material,  chemically  and  structurally Softwoods 
are  structurally  much  simpler  than  hardwoods  and  the  present  x^ork  is  mainly 
concerned  with  the  natural  multiphase  mixtures  occurring  in  Douglas  Fir 
and  Western  Hemlock.  For  purposes  of  analysis  a  simplified  diagrammatic 
sketch  of  a  softwood  sample  is  shown  in  Fig.  6.1 


Fig*  6.1  -  Diagrammatic  sketch  of  softwood  sample 


The  fiber  cavities  can  be  considered  to  be  ellipsoidal  in  shape  to  a 
good  approximation.  These  capillary  cavities  are  air  filled  up  to  the  fiber 
saturation  point  after  which  water  condenses  inside  the  cavities  starting 


. 


' 
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■with  the  smallest  capillary  diameter.  Below  fiber  saturation  the  water  is 
adsorbed  by  the  cell  walls  and  bonds  to  the  cellulose  fibers  in  layers  of 
one  to  seven  molecules  thick.  The  number  of  bonded  layers  increases  with 
increasing  moisture  content.  Stamm''  indicates  that  the  water  is  adsorbed 
on  the  surface  of  the  crystallites  which  are  cylindrical  in  shape  and  in 


Fig.  6.2  -  Specific  gravity  of  adsorbed  water  in  wood  vs.  moisture  content 


the  amorphous  regions  between  the  crystallites.  The  water  is  most  likely 

7 

adsorbed  as  little  spheres  ranging  in  curvature  from  .2  x  10  cm  to  20  x 

10"*^  cm.  Stamm  shows  a  curve  of  radius  of  curvature  of  the  adsorbed  water 

vs.  the  fractional  water  volume  in  the  cell  wall.  These  values  are  calcu- 

•><•**<* 

lated  from  Kelvin's  capillary  equation'  '  and  sorption  isotherms  for  the  wood 


*  Stamm  -  ref.  17 9  p.  48 
•vo*  Stamm  -  p.  390 


.  ' 


- 
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Table  6,1  -  Air  and  Water  Volumes  at  Different  Moisture  Contents 

and  Specific  Gravities 


%  Moisture  Content  %  WTater  Volume  %  Air  Volume  Specific  Gravity 


0 

0 

62.7 

5 

2.14 

61 . 4 

10 

4.32 

60.1 

15 

6.52 

58.7 

20 

8.69 

57.3 

25 

10,8 

56.0 

30 

12,8 

54.7 

0 

0 

64,9 

5 

2,01 

63.7 

10 

4.07 

62.4 

15 

6.15 

61.0 

20 

8.21 

59.7 

25 

10,2 

58.3 

30 

12,1 

57.1 

0 

0 

67.1 

5 

1.89 

65.9 

10 

3.83 

64.6 

15 

5.79 

63.3 

20 

7.74 

62.0 

25 

9.63 

60.7 

30 

11.4 

59.5 

0 

0 

69.3 

5 

1.76 

68.1 

10 

3.58 

66.9 

15 

5.43 

65.6 

20 

7.26 

64.3 

25 

9.05 

63.0 

30 

10.8 

61.8 

0 

0 

71.4 

5 

1.64 

70.3 

10 

3.34 

69.1 

15 

5.07 

67.9 

20 

6.79 

66 . 6 

25 

8.48 

65.4 

30 

10.1 

64.2 

0 

0 

73.5 

5 

1.53 

72.5 

10 

3.11 

71.3 

15 

4.72 

70.1 

20 

6.33 

68.9 

25 

7.9 

67.7 

30 

9.41 

66 . 6 

0 

0 

77.6 

5 

1.29 

76.7 

10 

2 . 64 

75.6 

15 

4.01 

74.6 

20 

5.4 

73.5 

25 

6.76 

72.4 

30 

8.06 

71.4 

.475 


.45 


.425 


.4 


.375 


.35 


.3 
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It  is  also  known  that  adsorbed  water  in  wood  is  compressed,  to  a  consider¬ 
able  extent.  The  specific  gravity  of  this  bound  water  ranges  from  1.3  at  0 

(22) 

moisture  content  to  1.115  at  the  fiber  saturation  point.  MacLean  gives 
a  curve  of  the  specific  gravity  of  adsorbed  water  vs.  moisture  content. 

A  least  squares  fit  to  the  data  given  by  MacLean  results  in  a  quadratic 
function  given  by,  (see  Fig.  6.2), 

fw  =  |.2S7-l.l04(Mq+  l.0966(MC)z  (6.1) 

where  ^  is  the  density  of  the  adsorbed  water  and  MC  is  the  percentage 

moisture  content.  RMS  deviation  from  the  given  data  was  less  than  .4%. 

(22)  (17) 

MacLean  and  Stamm  also  give  formulas  for  the  determination  of  the 
void  volume  and  the  specific  gravity  of  wood  as  functions  of  moisture  content. 
These  relations  are  repeated  here  for  convenience.  The  fractional  void 
volume  of  wood  can  be  calculated  from  equation  (6.2) 


V  =  1  —SI— _  4-  + 

va  1  +  loop  ^ 


1<46 


oop 


ik. 


(6.2) 


w  >°0 

where  is  the  moisture  content  above  fiber  saturation,  $>  is  the  density 
of  the  water,  5  is  the  specific  gravity  of  the  wood  at  a  given  moisture 
content,  based  on  ,  the  specific  gravity  of  wood  when  green.  The  equation 
for  S  is 

S  =  Sj[ I  -  .009 Sj  (30 -MC)] 

Finally  the  fractional  volume  of  water  is  given  as 


(6.3) 


Vw= 

100  f* 


(6.4) 


In  the  adsorption  range  of  wood,  roughly  0-30%  M;C.,  the  cell  walls 
swell  as  a  result  of  the  water  taken  up.  The  total  void  volume,  however, 


remains  relatively  constant,  see  table  6,1,  and  the  fiber  cavities  are 


filled  with  air  only. 


* 


Sizes  of  capillaries  and  crystallites  can  be  found  in  the  literature 
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• 

From  X-ray  diffraction  studies  the  crystalline  units  in  the  cellulose, 
referred  to  as  crystallites,  have  an  estimated  length  of  800  A°,  or  8  x  10  ^ 
cm,  and  an  average  cross  section  of  about  50  by  100  A° .  An  amorphous 

k 

sheath,  believed  to  exist  around  each  crystallite  ,  has  an  estimated  wall 
thickness  of  about  25  A° .  The  average  fiber  length  for  both  Douglas  Fir 

kk 

and  Western  Hemlock  is  3.15  mm.  Fibers  tend  to  be  rectangular  in  cross 

-4 

section  and  the  average  fiber  and  fiber  cavity  width  is  about  33.3  x  10  cm 

-4 

and  26.2  x  10  cm  respectively,  giving  a  fiber  wall  thickness  of  about 

— -j-  ^ 

7.1  x  10  cm  .  The  cell  walls  have  a  void  volume  less  than  2%.  These  are 
the  important  dimensions  necessary  to  predict  the  dielectric  properties  of 
the  wood  mixture  from  its  component  structure. 

On  basis  of  the  foregoing  discussion  it  is  reasonable  to  assume  that 
the  cell  walls  can  be  treated  as  a  homogeneous  solid  solution  of  cell  wall 
material  and  small  spheres  or  ellipsoids  of  adsorbed  water,  case  I.  Equally, 
one  can  postulate  that  the  crystallites,  impermeable  to  water,  are  surrounded 
by  an  almost  continuous  cylindrical  or  ellipsoidal  sheath  of  adsorbed  water, 
case  II.  Both  of  these  configurations  will  be  treated  in  the  following 
analysis . 

At  moisture  contents  above  the  fiber  saturation  point,  condensed  water 

bi 

in  the  capillary  structure  will  have  to  be  accounted  for.  There  is  evidence 
that  condensation  occurs  in  pit  membranes  on  the  cell  walls  at  relative  vapor 


* 

■k-k 

+ 

tt 


Stamm  -  p.  30 
Stamm  -  p.405 
Stamm  -  p.388 
Stamm  -  p.144 
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pressures  as  low  as  0.9  («=  20$  M.C.).  Initial  condensation  of  water  in  the 
pit  membranes  probably  consists  of  small  spherical  droplets.  Above  30$ 
the  water  condensed  in  fiber  cavities  will  assume  an  ellipsoidal  shape  and 
the  shape  of  the  air  inclusions  in  the  cavities  will  change  from  ellipsoidal 
to  circular  as  more  water  condenses  in  the  capillaries.  Unfortunately ,  it 
is  difficult  to  predict  how  much  air  is  held  as  spherical  inclusions  and 
how  much  is  held  in  ellipsoidal  inclusions. 

Case  I. 

A.  MIXTURE  DESCRIPTION 

It  is  assumed  that  the  soft  wood  cell  wall  material  is  a  mixture  of 
randomly  oriented  small  spherical  water  inclusions  dispersed  in  a  continuous 
medium  consisting  of  the  basic  wood  substance.  Within  this  host  medium, 
there  are  ordered  ellipsoidal  air  cavities.  Using  the  relations  given  in 
chapter  2,  it  is  possible  to  set  down  one  equation  describing  the  mixture 
dielectric  constant  of  the  3-phase  mixture  in  terms  of  its  component 
dielectric  constants  up  to  the  fiber  saturation  point  of  softwoods.  The 
ordered  nature  of  the  air  cavities  will  introduce  anisotropy  into  the 
mixture  as  a  whole.  Eventhough  the  basic  cellulose  structure  has  consi¬ 
derable  anisotropy  itself,  the  cell  wall  structure  is  such  that  the 
cellulose  fibers  must  be  considered  as  randomly  oriented  and  not  as 
ordered  units.  Hence,  the  anisotropy  of  the  mixture  due  to  the  cellulose 

( 5) 

fibers  is  assumed  negligible. 

The  total  air  volume  below  fiber  saturation  is  greater  than  50$  of 
the  mixture  volume  for  specific  gravities  of  .4  or  less.  However,  the 
interaction  effects  between  the  large  number  of  air  cavities  are  offset 
by  the  electrostatic  shielding  effect  of  the  high  dielectric  constant  of 


' 


, 
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the  cell-wall  mixture. 

i.  / 

B*  DIELECTRIC  PROPERTIES  OF  THE  MIXTURE  COMPONENTS. 

Wood  cell  substance  consists  basically  of  cellulose  and  Lignin , 

along  with  very  small  quantities  of  extractives.  Careful  measurements 
•  /0\ 

were  made  by  Trapp  and  Pungs'  '  on  the  dielectric  properties  of  cellulose 
and  Lignin,  as  well  as  on  the  air-free  cell  substance.  They  gave  the 
dielectric  constant  of  these  three  materials  as  a  function  of  frequency. 
The  low  and  high  frequency  dielectric  constant  as  well  as  the  relaxation 
time  of  the  cell  substance  are  estimated  from  their  data  as  6.6,  2,  and 
2  x  10^  sec  respectively  at  20°C.  Similarly,  for  water  the  low  and 
high  frequency  dielectric  constants  are  given  by  Hasted (51)  as  80.36, 

5.5,  with  a  relaxation  time  of  .923  x  10^*  sec  at  20  C. 

However,  when  water  is  molecular ly  bound  to  the  host  material,  as 

is  the  case  for  wood  below  fiber  saturation,  the  dielectric  properties 

of  the  water  can  be  considerably  different  from  the  pure  liquid  water 

( 59) 

values.  Haggis,  Hasted,  and  Buchanan  showed  this  quite  clearly  for 

water  in  solutions.  Windle  and  Shaw^26;  actually  calculated  a  value  for 

the  dielectric  constant  of  localized  and  mobile''  water  in  wool-water 

systems  from  their  experimental  data.  The  latter  authors  assume  two 

states  of  adsorbed  water  whereas  experimental  evidence  presented  in  this 

thesis,  see  for  example  Fig.  5.26,  and  arguments  given  by  McLaren  and 

Rowen^65^  in  their  revie?/  paper  on  water  sorption  shows  that  mobile 

¥  "Ijodal  17, pel "water  refers  to  that  water  which  is  attached  directly  to 
polar  groups  in  the  host  material  giving  rise  to  high  activation 
energies  of  polarization*  Mobile  water  refers  to  thau  water  which  is 
indirectly  attached  to  the  mixture  through  other  molecules,  hence, 
having  lower  activation  energies  of  polarization. 
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water  likely  has  dielectric  properties  varying  continuously  between 
the  properties  of  localised  water  and  liquid  wa ter* 


The  dielectric  constants  of  bound  water  are  not  adequately  known 
nor  readily  calculated.  By  combining  the  experimental  data  on  Douglas 
Fir  as  given  in  chapter  5  and  the  mixture  relations  developed  in 
chapter  2,  it  is  possible  to  determine  semi-empirical  estimates  for  the 
dielectric  constant  of  adsorbed  water  as  a  function  of  moisture  content. 


C.  ACTUAL  IvUXTURE  RELATIONS  USED  FOR  CASE  I. 

Under  the  assumptions  of  section  A  a  3-phase  mixture  relation  is 

to  be  found  for  ordered  ellipsoidal  air  cavities  and  random  spheres  of 

water  dispersed  in  a  homogeneous  cell  wall  material  as  shown  in  Fig.  6.3. 

The  internal  field  for  the  water  spheres  can  be  found  from  equation  (2.66) 

1  * 

averaged  over  all  orientations  using  n=-  with  6,  ='£w  ,  =  e,v  gives 


■y 


c 


tv 


d 


w 


(6.5) 


where  £wis  the  dielectric  constant  of  the  adsorbed  water  and  is  the 
effective  dielectric  constant  of  the  material  surrounding  the  water 
inclusions.  For  the  air  inclusions  equation  (2.66)  is  used  again  with 
£,  =  £*,  £1  =  £*  which  yields 


c  *  r  p* 


+  nja.-d)] 


(6.6) 


where  £0  is  the  dielectric  constant  of  air,  £*  is  the  effective 
dielectric  constant  of  the  material  surrounding  the  air  inclusions. 

Substitution  of  equations  (6.5)  and  (6.6)  into  (2.22),  using  6^  =  P c. 
gives  the  3-phase  mixture  relation 
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— £c)£q +  3VW  few  -kU* 

£<*  +  na  (£„-£*)  2-C+£„ 


(6.7) 


vAere  £mis  the  mixture  dielectric  constant,  £c  is  the  dielectric 
constant-  of  the  dry  cell  wall  material,  Va  and  are  the  fractional 
volume  filling  factors  of  the  air  and  water  respectively. 


Ellipsoidal  ait  cavities  -  £ 


a. 


0*"/ Spherical  water  inclusions  -  £ 
Cellwall  material  -  £c 


Fig*  6*3  -  Cellulose-air-water  mixture. 


The  maximum  and  minimum  values  that  and  £a  can  have  lie 
between  the  maximum  and  minimum  values  of  £m  and  £c  *  In  the  following 
results  it  has  been  assumed  that  =£a-£w  since  the  lower  limit  given  by 
£c  gives  dielectric  constants  far  below  the  experimental  values.  From 
the  fiber  cavity  dimensions  for  softwoods  it  is  found  that  in  the  direction 
of  the  major  axis  of  the  air  inclusions  n-0and  in  the  other  two  directions 
n=-J-,  The  standard  formula  for  depolarization  coefficients  was  used  to 
obtain  these  values  (see  equation  (2,50)  and  also  reference  2),  Assuming 
a  specific  gravity  of  .400,  VA  and  V;v  can  be  found  from  Table  6.1,  The 
dielectric  constant  of  the  cellwall  material  is  taken  from  Trapp  and 
Pungs7data  as  £c=3.5  -  j.012  at  2450  MHz*  Different  values  for  £w  are 
used  in  the  following  results  to  indicate  the  fact  that  the  adsorbed 
water  behaves  different  than  pure  liquid  water.  Plots  of  the  dielectric 
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constant  and  dielectric  loss  versus  moisture  content  at  20  are  given 
in  figures  6*4  and  6*5  under  the  assumption  of  spherical  water  inclusions* 
Experimental  data,  taken  from  figures  5*13  and  5.14  are  given  for  com¬ 
parison  and  are  tabulated  for  convenience  in  Table  6.2# 


TABLE  6*2  -  Experimental  values  of  the  dielectric  constant  of  long¬ 
itudinal  Douglas  Fir  at  20  C  and  2450  MHz  (s.g.  *  .4). 


Moisture 
Content,  % 


0 

1.78 

.04 

5 

2.14 

.14 

10 

2.5 

.325 

15 

2.9 

.58 

20 

3.6 

.86 

25 

4.35 

1.16 

30 

5.15 

1.46 

The  parameter  in  figures  6.4  and  6.5  is  the  dielectric  constant 
of  the  included  water.  Using  the  value  for  pure  liquid  water gives 
a  complex  dielectric  constant  which  falls  far  below  the  experimental 
data.  Similarly  the  value  for  mobile  water  calculated  by  Windle  and 
Shaw^*^ ,  namely  £w=60-j  12.5  gives  values  much  lower  than  the  exper¬ 
imental.  Even  values  of  the  dielectric  constant  of  adsorbed  water  as 
determined  semi-empiric ally  in  the  following  section  give  values  which 
are  low.  The  only  remaining  cause  for  the  discrepancy  is  that  the 
assumed  model  for  the  mixture  is  wrong.  Since  the  air  inclusions  are 
known  to  be  ellipsoidal  in  shape  from  microscope  studies,  the  assumed 
shape  of  the  water  inclusions  must  be  the  cause  of  the  disagreement. 

It  can  further  be  stated  that  the  shape  of  the  water  inclusions  must 
be  such  as  to  give  rise  to  a  very  small  reaction  or  depolarization 
field  to  cause  the  required  increase  in  dielectric  constant,  for  ex¬ 
ample  f ellipsoids  or  ellipsoidal  shells. 


■  . 
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Fig.  6.4  -  Theoretical  dielectric  constant  of  a  softwood 

mixture  versus  moisture  content  with  the  dielectric 
constant  of  included  water  as  a  parameter;  E-field 
in  the  longitudinal  grain  direction  at  2450  MHz  for 
a  s.g.  of  .4;  spherical  water  inclusions. 
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Fig.  6.5  -  Theoretical  dielectric  loss  of  a  softwood  mixture 

versus  moisture  content  with  the  dielectric  constant 
of  adsorbed  water  as  a  parameter;  E-field  in  the 
longitudinal  grain  direction  at  2450  MHz  for  a 
s.g.  of  .4;  spherical  water  inclusions. 
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Case  II. 

A.  MIXTURE  DESCRIPTION 

In  this  case  it  is  assumed  that  the  included  water  occurs  in  the 
form  of  ellipsoidal  shells  and  the  mixture  further  contains  ordered  air 
cavities  as  in  Case  I.  The  arguments  regarding  anisotropy,  localized 
and  mobile  water  for  Case  I  apply  equally  well  here  and  will  not  be 
repeated.  The  mixture  model  used  here  is  depicted  in  figure  6.6. 


Ordered  Ellipsoidal  air 
cavities  -  £a 

Random  Ellipsoidal  water  - 
shells 

Cellwall  material  -  £c 


£|v 


Fig.  6.6  -  Cellulose-air-water  mixture;  Case  II. 


B.  MIXTURE  RELATIONS  USED  IN  CASE  II. 

In  a  manner  very  similar  to  that  given  under  Case  I,  the  appropriate 
internal  fields  are  given  by  equations  (2.65)  and  (2.66)  for  the  water 


and  air  inclusions  respectively .  Substitution  of  these  internal  field 
expressions  for  Tt  in  equation  (2.22)  then  gives  the  final  equation  for 
this  3-phase  mixture. 


(6.8) 
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■where  and  are  the  dielectric  constants  of  air  and  of  the 
material  surrounding  the  air  inclusions ,  Ho,  is  the  depolarization 
coefficient  for  the  air  cavities  which  is  vjhen  the  E-field  is 

applied  in  the  direction  of  the  major  axis.  The  second  term  concerns 
the  random  water  inclusions.  It  is  assumed  that  the  material  in 
region  three,  the  inner  ellipsoid  is  that  of  the  dry  cell  wall  (£j  =  £c), 
the  material  in  region  1  is  assumed  to  have  a  dielectric  constant 
equal  to  that  of  the  complete  mixture,  i.e.,  £,  =  £wand  £2  is  the 
dielectric  constant  of  the  adsorbed  water  (£2.-^)  •  The  host  material 
has  a  dielectric  constant  equal  to  that  of  the  dry  cell  wall  (£y,=  £c) 
and  the  depolarization  coefficients  Y\w  and  are  defined  by 
equations  (2.50)  and  (2.52)  respectively.  These  coeeficients  can  be 
calculated  using  the  dimensions  of  the  shell  structure. 

As  indicated  in  the  beginning  of  this  chapter,  little  is  known 
about  the  exact  sizes  of  the  amorphous  and  crystallite  regions. 

Hw  and  are  based  on  the  approximate  amorphous  sheath  dimensions 
given  earlier  in  this  chapter,  namely,  the  axes  a,  b  and  c  of  the 
out e  ^ellipsoid  are  approximately,  425,  50,  75  A  and  for  the  inner 

Vf 

ellipsoid  a* ,  b»  and  c*  are  400,  25,  and  50  A°  respectively. 

These  values  were  substituted  into  equation  (6.8)  which  was 
then  solved  numerically  for  two  different  values  of  the  dielectric 
constant  of  the  adsorbed  water.  The  results  are  presented  in 
figure  6.7.  Experimental  values  for  the  dielectric  constant  of 

*  It  turns  out  that  the  depolarization  coefficients  along  the  major 
axis  are  nearly  equal  to  zero  and  along  the  two  minor  axes  are  ap¬ 
proximately  equal  to 
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Fig.  6.7  -  Complex  dielectric  constant  of  a  cellulose- air-water 
mixture  versus  moisture  content ,  assuming  ellipsoidal 
shell  water  inclusions  with  experimental  data  for 
longitudinal  Douglas  Fir  given  for  comparison; 

{=  2450  MHz,  s.g .=  .400,  T=20°C. 


Dielectric  loss  factor  £  n 


- 


140 


longitudinal  Douglas  Fir  are  given  for  comparison.  It  is  observed 
that  the  calculated  dielectric  loss  is  too  low  and  the  real  part  of 
the  dielectric  constant  is  too  high  at  low  moisture  contents  compared 
to  the  experimental  data. 

Ope  probable  reason  for  the  disagreement  with  the  mixture  models 
given  so  far  is  that  the  assumed  dielectric  constant  of  the  adsorbed 
vrater  is  incorrect.  To  resolve  this  question  an  estimate  of  the  bound 
water  dielectric  constant  was  obtained  by  the  following  procedure. 

The  wood  is  assumed  to  consist  of  a  homogeneous  host  material  of 

■55- 

dielectric  constant  in  which  ellipsoidal  air  cavities  are  included. 
Since  the  depolarization  coefficient  for  the  air  cavities  is  nearly 
equal  to  zero  along  the  major  axis,  it  is  found  that  equation  (2.66) 

substituted  into  (2.22)  simplifies  to 
—  £-m,  t  ( £*  "  £> n,) 

Hence,  the  dielectric  constant  of  the  ‘cellulose-water *  host  material 
is  given  by 

£m,=  (6*9) 
i  -va  ' 

A  plot  of  equation  (6.9)  is  given  in  fig.  6.8.  The  dielectric 
constant  £.  thus  obtained  is  the  value  for  a  mixture  of  water  inclu- 

m, 

sions  in  the  basic  cell  wall  material.  The  value  of  £c  the  cell  wall 
material  is  about  3.5  -j.012  and  the  dielectric  constant  of  the 
adsorbed  water,  can  be  assumed  to  be  considerably  greater  than  the  value 

*  Since  the  water  inclusions  are  orders  of  magnitude  smaller  than  the 
air  cavities  below  fiber  saturation,  this  assumption  is  valid. 
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Fig.  6.8  -  Dielectric  constant  of  the  cell  wall  water  mixture 
with  no  air  inclusions,  versus  moisture  content, 
calculated  from  equation  (6.9)  using  the  experi¬ 
mental  dielectric  constant  of  Douglas  Fir  from 
Table  6.2;  -f=  2450  MHz,  s.g.=  .400,  T=20°C. 
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ol  £c  .  The  final  results  will  show  this  assumption  to  be  valid, 
hence,  £c  ^  Lw .  To  simplify  the  calculation,  the  water  is  assumed 
to  oe  in  the  form  of  thin  ellipsoidal  shells  so  that  Y\w—Y\ J  and 
V=V  .  Under  these  conditions,  equation  (6.8),  remembering  VA  =0, 
and  is  now  >  simplifies  and  when  solved  for  the  dielectric 
constant  of  the  adsorbed  water  gives 

~  ~  £c  ( 1  ~~Vw)  (6.10) 

Vw 

where  £m>  is  given  by  equation  (6.9). 

The  numerical  results  for  this  approximation  to  the  bound  water 
dielectric  constant  are  presented  graphically  in  figure  6.9.  It  is 
clear  that  the  adsorbed  water  behaves  differently  than  pure  water. 
Windie  and  Shaw’s  value  for  £w,  60- j  12.5  at  3  GHz  seems  to  agree 
reasonably  with  the  values  calculated  here.  However,  the  dielectric 
loss  of  the  adsorbed  water  as  given  here  is  a  lot  higher  than  theirs. 
If  now  the  values  of  £w  from  figure  6.9  are  used  in  equation  (6.8) 
to  calculate  the  total  mixture  dielectric  constant,  assuming  the 
amorphous  sheath  model  again,  it  is  found  that  agreement  between  the 
experimental  and  theoretical  values  is  very  good,  indicating  that 
the  approximations  made  in  deriving  equation  6.10  were  reasonable. 
Figure  6.10  shows  the  resultant  theoretical  and  experimental  curves 
as  a  function  of  moisture  content. 

The  calculation  is  repeated  with  *V=J-  and  the  result  plotted  in 

*  This  may  seem  to  be  a  crude  approximation,  but  it  is  shown  to  give 
reasonable  results  and  the  errors  due  to  the  approximation  can  be 
shown  to  be  about  25$. 
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%  Moisture  Content 


-  Approximate  value  for  the  dielectric  constant  of 
adsorbed  water  versus  moisture  content  based  on 
equation  (6.10)  assuming  random  thin  shell  inclu¬ 
sions;  f  =  2450  MHz,  T=20°C. 


Fig.  6.9 
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Fig.  6.10  -  Theoretical  and  experimental  plot  of  the  dielectric 
constant  of  Douglas  Fir  versus  moisture  content 
with  the  E-field  in  the  longitudinal  grain  direc¬ 
tion  using  semi-empirical  values  for  the  adsorbed 
water  dielectric  constant;  Mixture  model:  ellip¬ 
soidal  water  shells  and  ellipsoidal  air  cavities; 
f=  2450  MHz,  T=20°C.,  s.g.=.400. 
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Fig.  6.11  -  Theoretical  and  experimental  dielectric  constant 

of  Douglas  Fir  versus  moisture  content  mth  E-field 
in  the  tangential  grain  direction;  f  =  2450  MHz, 
s.g.«.400,  T=20°C. 
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figure  6,11.  This  plot  thus  gives  the  dielectric  constant  in  the 
radial  or  tangential  direction  since  the  depolarization  coefficients 
were  assumed  equal  along  the  two  minor  axes.  For  comparison  the 
experimental  values  for  tangential  Douglas  Fir ^  given  in  chapter  5? 
are  superimposed  and  the  agreement  is  fair.  Therefore 5  the  assumption 
that  the  anisotropy  of  'wood  is  mainly  due  to  the  ordered  nature  ox 
the  air  cavities  is  shown  to  be  reasonable  also. 


Finally  the  static  and  high  frequency  dielectric  constants  can 
be  calculated  by  substituting  the  corresponding  low  and  high  frequency 
component  dielectric  constants  (found  under  case  I  section  B  in  this 
chapter)  into  equation  (6.0),  Figure  6,12.  shows  that  hero  also  reas¬ 
onable  agreement  between  theoretical  and  experimental  values  is  found. 


It  can  thus  be  concluded  that  the  cellulose-air-water  mixture 
as  found  in  wood  is  adequately  described  by  the  amorphous  sheath 
approximation.  Since  the  sheath's  dimensions  are  only  approximately 
known  the  results  could  be  improved  by  a  better  understanding  and 
description  of  the  actual  crystalline  -  amorphous  structure  of  the 
cell-avail  material.  The  estimate  of  the  bound  water  dielectric 
constant  found  here  should  be  a  useful  tool  in  the  analysis  of  other 
structures  containing  adsorbed  water.  More  detailed  description  cf 
the  geometry  is  essential  for  further  improvement  in  the  values  given 
here , 


Thus  the  complicated  wood  mixture  analyzed  in  this  thesis  can 
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Fig.  6.12  -  Theoretical  and  experimental  lovr  and  high  frequency- 
dielectric  constants  versus  moisture  content  for 
Douglas  Fir  -with  E-field  in  the  longitudinal  grain 
direction;  f =  2450  MHz,  s.g#=.400,  T=20°C. 


adequately  be  described  by  the  new  mixture  relation  using  the  shell 
inclusions  and  the  semi-empirical  values  of  the  bound  water  dielec tri 
constant  as  calculated  in  this  chapter 0  Previous  mixture  relations 
were  unable  to  describe  the  wood  mixture  satisfactorily  as  a  function 


of  moisture  content. 


CHAPTER  7 


CONCLUSION 


A  multiphase  dielectric  mixture  theory  is  presented  vihich  allow*:; 
analysis  of  ellipsoidal  shell  inclusions.  The  mixture  equations  developed 
by  others  for  ellipsoidal  and  spherical  inclusions  are  shown  to  be  a 
subset  of  the  more  general  solution  presented  in  this  thesis  *  Examples 
of  shell-type  inclusions  can  be  found  in  natural  materials  such  as 
brick  and  wood .  The  latter  case  is  analyzed  by  moans  of  the  mixture 
theory  given  here  and  the  results  indicate  that  the  model  chosen  for 
the  mixture  is  reasonable  *  Dielectric  properties  for  adsorbed  rater 
are  determined  semi-empiric  ally  from  experimental  Douglas  Fir  data 
which  indicate  that  adsorbed  water  behaves  very  differently  than  purs 
liquid  water .  The  numerical  values  given  here  for  the  bound  voter 


e  stir; 

latcs 

since 

the  results 

depend 

chose:' 

<  for 

-.rood , 

Tills  model 

is  not 

necessarily  exact  because  relatively  little  is  known  abort  the  way 
water  is  held  in  the  cell  walls „ 


Behavior  of  mixed  dielectrics  as  a  function  of  temperature  and 
moisture  content  is  analyzed  under  Inc  assumption  In  at  disi/rwuit—ons 


ouglas  Fir  and  Vie  stem  Hemlock  confirm  the  theory  to  better 
f  As  a  result  it  is  possible  to  describe  the  dielectric 


and  frequency  knowing  only  the  activation  energy  of  polarization  and 


moi 

sture  i 

cf 

relax  a 

dat 

a  on  B 

tha 

1  •>’ 
j _ i 

b 

viropertie 

150 


the  static  and  high  frequency  dielectric  constant  of  the  material  as 
a  function  of  moisture  content.  If  easier  methods  of  measurement  could  be 
found  to  determine  these  three  parameters ,  it  would  greatly  facilitate 
the  acquisition  of  dielectric  data  yet  needed  for  many  materials. 

It  is  evident  from  the  experimental  method  chosen  for  this  work  and  the 
material  under  test  that  acquiring  enough  data  to  adequately  describe 
the  dielectric  properties  is  a  laborious  and  time  consuming  process. 
Nevertheless,  the  data  obtained  in  this  manner  allow  a  good  description 
of  a  rather  Yri.de  class  of  cellulose  containing  materials.  The  basic 
accuracy  of  the  measurement  technique  used  is  extremely  good,  but  if 
all  possible  variations  due  to  the  material,  control  and  readout 
devices  used  are  taken  into  account  the  overall  accuracy  of  the 
measurement  deteriorates  to  +  10$. 

Methods  described  and  data  obtained  in  this  work  can  be  used  to 
advantage  in  the  description  of  microwave  heating  and  drying  processes. 
From  the  results  on  wood  it  is  possible  to  select  optimum  operating 
frequencies  and  temperatures  for  such  processes.  Applying  the  results 
to  control  instrumentation  design  allows  one  to  select  a  frequency 
which  minimizes  calibration  problems,  by  viiich  present  instruments, 
for  example,  moisture  meters  are  plagued.  To  adequately  design 
processing  systems  and  accompanying  instrumentation  requires  a 
thorough  knoYriLedge  of  the  materials  to  be  treated. 


This  v/ork  has  laid  the  basis  for  a  better  understanding  of  the 
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properties  of  adsorbed  water  and  the  way  in  which  it  occurs  in 
cellulose  materials.  Electron  microscope  and  X-ray  diffraction 
studies  coupled  with  more  detailed  dielectric  measurements  will 
likely  resolve  the  remaining  uncertainties  in  the  model  of  cel- 
lulos e-air-water  mixtures  given  here.  Complexity  of  the  result¬ 
ing  mathematics  is  greatly  offset  by  modem  computing  facilities 
which  extend  the  research  power  of  a  scientist  by  orders  of  mag¬ 
nitude  , 
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APPENDIX  A 


.U  rv 


ERROR  .ANALYSIS 

In  this  section,  the  basic  equations  determining  the.  errors  in  the 
dielectric  constant  are  worked  out  and  some  useful  approximations  are  stated. 
The  results  are  applied  to  the  dielectric  measurement  techniques  used  in 
Chapter  A. 

Let  the  dependent  variable,  E,  and  the  independent  variable,  V,  be 
related  by 


n 


(V) 


A  ( 1  ') 


E  and  V  are  any  two  variables  and  have  no  other  significance.  Let  £E  be 
the  increment  in  E  due  to  an  incremental  change, A  V,  in  V.  Then  by 
definition  of  the  derivative, 


Urn  aK/M'=  d  E/d  V 

A  V*  0 


(2) 


When  A  V  is  very  small 

At  =  (d B/dV)  AV 


A  (3) 


However ,  the  exact  expression  f or  A E  can  be  obtained  by  a  Taylor  series 


expansion  of  equation  X2 j 


y  (36) 


Thu 


AF 


f/’./j  /  \/ } 
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A  (4) 


hence , 
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AE~  { '( Vj  A  V  +  f}(V)(AV)l/z!  + 

,  .  AK->) 

+  f3iy)(A\i]/3!  +  •  •  • 

J  ■> 

Where  •f(V) }  rf  V,'  ,  etc.,  designate  the  first,  second  and  higher  order 

derivatives  of  f  w.r.t.  V.  More  generally,  E  can  be  a  function  of  more 
than  one  variable,  in  other  words,  V  can  be  a  vector,  VI,  V2,  ... 
Similarly,  E  itself  can  be  a  vector  of  dependent  variables,  El,  E2,  ..., 
in  which  case  the  functional  relationship  can  be  written  as, 


A  (6) 


*> 

*• 


Writing  the  incremental  changes  in  E  due  to  incremental  changes  in  V  by 
using  a  Taylor  series  expansion  for  n  independent  variables,  equation  A (6) 


becomes 
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Where 


^  i 


Yi 


and  so  on . 


If  AVI,  A V2  ,  •  •  •  are  small,  as  is  usually  the  case,  higher  order  terms 
can  be  neglected  and  hence  A (7)  can  be  rewritten  as, 
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Equation  £(7),  therefore,  gives  the  total  change  in  the  dependent  variables, 
Ei,  E2,  as  a  function  of  the  change  in  the  independent  variables,  VI,  V2,... 

In  practice,  equation  A  (8)  becomes  a  very  good  approximation,  and  usually 
consists  of  just,  a  few  terms,  since  many  of  the  derivatives  will  be  zero 
and  the  increments  are  usually  small , 

Consider  now  a  particular  measurement  technique,  for  example,  the 
bridge  method  used  in  chapter  -•  for  measuring  dielectric  constants.  The 
dependent  variables  are  the  real  and  imaginary  parts  of  the  dielectric 
constant,  £  '  and  6  M  respectively.  For  the  case  of  guided,  fundamental  mode, 
electromagnetic  propagation,  the  functional  relationships  are  given  by- 
equations  A (9)  and  A (10) 


£ 1  ~  p  ■+  ( )\0  fi  n )  ( s  c  ‘ ) 

£“=2(\0/27  rFocP  .  A(10) 

where  ©i  and  j3  ,  the  two  independent  variables,  are  the  attenuation  and  phase 
constants  of  the  material,  ac, is  the  f r ee-spac e  wavelength,  and  p  is  the 
ratio  squared  of  the  cutoff  to  signal  frequency.  By  application  of  equation 
A(7)  it  is  found  that, 


AC  ID 
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since  all  other  derivatives  vanish.  Evaluating  the  partials  by  using  equa¬ 
tions  A(9)  and  A(10)  and  neglecting  second  order  terms  one  obtains 

AV~  2(A/2T[f  (B/j!? 


A(X2) 
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Thus,  equation  A, (12)  gives  the  error  in  dielectric  constant  directly  in 
terms  of  the  two  variables  of  and  B  which  in  turn  are  to  be  determined 
exp  er  ini  en  tally. 

Note  that  c<  and  B  are  themselves  functions  of  some  signal  amplitude, 
r,  and  phase  angle,©  .  For  the  bridge  method  under  discussion  here,  these 
functional  relationships  are  known  and  can  be  written  as,  sec  chapter  4, 

of-  (I  /2  s)  cosh  ’{ (r2/4)  +  Vll-r1//^  +  r2  s  in2  e  |  A'- 


( i /z -()  cos'1 1  (rzM)~\/[i-ritt):L  +  r2sin29l 


A(14) 


where  is  the  sample  length,  F  is  the  amplitude  ratio  of  the  detected 
signal  at  two  different  positions,  and  ©  is  the  phase  angle  between  the 
two  signals.  By  again  using  equation  A(7)  we  can  obtain  the  incremental 
change  in  X  and  B  due  to  small  changes  in  F  and  £  .  Hence, 

B/a  r  )t  r + {'<>  f  /o 9i l  o  4(1/2)  |  (>>'■  e/e r%  rf + 
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+  2(oV/erS5)/-rA0  (av./302)(2i0;~J- 

Again  second  and  higher  order  terms  can  be  neglected  since  4!'  and  L  9  are 
small  if  a  reasonably  good  attenuator  and  phase  shifter  are  used  in  the 
experiment . 


Evaluating  .dv  5r  from  equation  A(13)  gives 


,  1 
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It  can  be  shown ,  remembering  that  ;  is  less  than  unity,  that. 
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and  similarly, 


0  B  jb  y  -  5/2  ■£ 

do </dQ  £  r/2  i 


A  (19) 
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Using  these  upper  error  limits  in  equations  A  (15)  and  A (16)  yield: 

/i  B  o) 

/idiC  (i  \l)  [a  r  +  (f/?){i  el 

If  the  experimental  errors  A Y  and  A  0  are  known,  then  the  error  in  and  £  " 
can  be  computed  by  use  of  equation  A(22)  and  A(13)  for  any  value  of c<  and  ;  . 

As  an  example,  some  typical  values  will  be  substituted  in  equation  A(12) 
The  resulting  errors  are  indicative  of  the  overall,  accuracy  of  this  measure- 
men  t  t echniqu e . 


The  uncertainty  in  the  attenuation  measurement  is  ±  . 05db  and  the  phase 


shift  error  is  ±  2  degrees.  Further  for  a  loss  tangent  of  .  ,  :  6  .5  for  a 
half  wavelength  sample.  For  X  =12 cm,  :  AACli),  2cr.  .  By  definition, 

r -Emin/ Umax ,  and  A-20  log  .r.  Thus,  by  equation  A(3) 

X  vJ 

/■  r~  frj8,G85)Ar'\  -1.0021 
AO  ~  2/3 SO  -  1  .OOS'6 

Using  equation  A(22)  the  errors  in  B  and  ■  are 

AB£  ±  .002 .4 
A  &  A  +  .00. M3 

Finally  using  equation  A(12)  the  errors,  due  to  variations  ?.n  r  and  o  ,  at 
fixed  temperature,  moisture  content,  etc.,  are 

A  £’5-  1 .02?. 

A  +  .0123 

Yov  a  dielectric  constant  of  2  and  loss  tangent  of  1  the  percentage  ct:o?.s 
in  f  ■>  and  l''  are  less  than  ±1%.  Similar  calculations  for  a  dielectric 


' 


*  i 


constant  of  2  and  a  loss  tangent  of  .01  yields  percentage 
£ !t  which  are  less  than  ±  0.1%. 

For  an  error  of  .010"  or  .025  cm.  in  the  measurement  of  t? 
errors  in  £ ’  and  £ "  are  about  ±1%,  From  experimental  dat 
ship  between  o'  and  temperature,  T,  was  obtained  and  thus 
and  £  "  w.r.t.  temperature  were  calculated.  It  was  four: 
to  temperature  variation  were  always  less  than  ±1%. 
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errors  in  £ ’  and 

the  corresponding 
a  the  relation- 
the  errors  in  £ ’ 
d  that  errors  due 
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